THE  DEVELOPMENT  OF  A COMPUTER  SPEECH  PROCESSING  SYSTEM  AND  ITS  USE  FOR 
THE  STUDY  AND  DEVELOPMENT  OF  PROCESSING  METHODS  FOR  ENHANCING  THE 
INTELLIGIBILITY  OF  SPEECH  IN  NOISE 

Russell  J.  Niederjohn 
Marquette  University 

Robert  A.  Curtis 

Rome  Air  Development  Center 


G8ff  WMLWtt  TO  1)06 
CTT  FULLY  LEGIBLE  PRGi&il&H 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whit  D*f  Enfrod) 


REPORT  DOCUMENTATION  PAGE 


2,  GOVT  ACCESSION  NO. 


THEJJEVELOPMENT  OF  A ^COMPUTER  SPEECH  PROCESSING 
.SYSTEM  AND  ITSJJSE  FOR  THE, (STUDY  AND  DEVELOPMENT 
OF.^JROCESSING, METHODS  FOR  ENHANCING  THE, INTELLI- 


GIBILITY OF,  SPEECH  IN,  NOISE  • 

AUINURIW 


Russell  Jy^Niederjohn  Robert 

TBiquULLt!  fflltvuigity  ■ (KADC7T 
Milwaukee,  Wisconsin  53233  Griffiss  AFB  NY  13441 


A^Curtls  ^ 


» PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Rome  Air  Development  Center  (IRAD) 
Griffiss  AFB  NY  13441 


II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Same 


U monitoring  agency  nau^kXqobls 
Same  l / ✓ 


S (II  ditto rut  from  Controlling  O ttlco) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIEN 


aeJint' 


ATALOG  NUMBER 


S.  TVPF  n rypfiBT  a PFBinn 

In-houseJ Interim  jfep'Skt 
6 Jun  0-26  Aug  77; 

«.  FLHIUWMWU  URU.  RERUN  I 


N/A 


6.  CONTRACT  OR  GRANT  NUMBER!*; 


N/A 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A WORK  UNIT  NUMBERS 


310U.G 

70550707 


IS.  NUMBER  OF  PAGES 


■lM- 

re.  S( 


SECURITY  CLASS,  (ol  l hit  report) 

UNCLASSIFIED 


IS*.  OECL  ASSI  FI  CATION/ OOWNGR  AOING 
SCHEDULE 

N/A 


16.  DISTRIBUTION  STATEMENT  (ol  (III*  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  Ml*  ebetrect  entered  In  Bloch  20,  II  dltlerent  It  mi  Report) 

Same 


16.  supplementary  notes 

RADC  Project  Engineer:  Captain  Robert  A.  Curtis  (IRAD) 


19.  KEY  WOROS  ( Contlnu*  on  reverse  oldo  It  necotoory  end  Idenf/fy  by  block  number) 

Speech  Enhancement 
Speech  Intelligibility 
Computer  Speech  Processing 


20.  ABSTRACT  (Continue  on  reveree  e/de  It  neeeesery  end  idontily  by  block  numbor) 

The  enhancement  of  speech  intelligibility  in  the  presence  of  high  noise 
levels  is  an  importantjeonsideration  in  many  practical  communication  situations. 
The  work  documented  in^this  report  describes  'the^ results  of  an  examination  of 
mvmsV  methods  for  processing  speech  so  as  to  enhance  its  intelligibility  in 
r%1ie  ^ resente  of  wideband  random  noise"^  ,/our  par44cular- pro cessing  methods  were: 
considered:  faj-y  INTEL,  a method  which  involves  processing  in  both  the  first  and 
second  order  spectral  domains;  (b)  spectral  subtraction,  a method^ which  involves 
(over)  -lC. 


DD  /jS™,  1473  EOITION  OF  I NOV  61  IS  OBSOLETE 


UNCLASSIFIED 


^0  9 JfO 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whon  Do to  Entorod) 


C 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGEfW7i»n  Q«l«  Enfnd) 

/ . J £ / A ' 3 


/ > a simple  subtraction  of  the  average  noise  spectrum  from  the  first-order  spectrum; 

jninimum  mean  square  error  filtering,  <a. -aathod\hich  involves  filtering  speech 
in  such  a way  as  to  minimize  the  mean  square  error  between  a signal  and  its 
expected  value  in  noise;  andjfi)  methods  based  upon  suppressing  the  frequency 
content  of  a speech  plus  noise  signal  between  pitch  harmonics  of  the  speech 
signal. 

^>To  carry  out  a study  of  methods  to  enhance  speech  intelligibility  in  noise, 
two^ general-purpose  computer  processing  systems  were  implemented.  The  first  ? 
'vas',  a terminal  interactive  system  for  'the?  generation,  analysis,  and  graphic 
display  of  synthetic  voiced  speech  sounds,  -,Through  the^'use  of  this  system,  a 
> considerable  insight  into  the  effect  of  various  processing  algorithms  upon 
speech  and  upon  speech  in  noise, has  been  effected. 

■^The  second  computer  processing  systMi  has  been  developed  for  the  processing 
of  real  speech. J^Thia  system* involves  the- use  of  a DDP-116  data  converter  and  a 
Honeywell  6000  Computer .^Communication  between  these  two  computers  is  by  means 
of  seven  track  magnetic  tape.  In  use,  this  system  facilitates  the  input,  pro- 
cess, and  playback  of  real  speech  utterances.  Through  this  system,  the  effect 
of  numerous  processing  algorithms  upon  normal  speech  in  noise  has  been  studied. 

While  both  of  these  processing  systems  have  been  developed  for,  and  applied 
to,  a study  of  processing  techniques  for  enhancing  the  intelligibility  of  speech 
in  noise,  the  computer  programs  generated  have  purposely  been  made  general  pur- 
pose so  as  to  facilitate  their  future  use  at  RADC  for  other  speech  processing 
and  signal  processing  tasks. 

The  effect  of  several  processing  algorithms  (based  upon  the  four  methods 
mentioned)  has  been  studied  for  numerous  synthetic  voiced  speech  sounds  and  for 
two  12-second  real  speech  utterances.  These  two  speech  utterances  were  generated 
by  a male  talker:  one  utterance  in  a signal-to-noise  ratio  of  +6  dB,  the  other 
utterance  in  a signal-to-noise  ratio  of  -6  dB.  Overall  results  indicate  that 
while  the  greatest  speech  enhancement  success  has  been  achieved  with  the  INTEL 
and  the  minimum  mean  square  error  filtering  methods,  the  four  methods  studied 
each  offer  a significant  potential  for  speech  intelligibility  enhancement  in 
noise. 
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Figure  1 - An  illustration  of  the  DDP-116/Honeywell  6000  magnetic- 
tape  communication  method. 

Figure  2 - An  overview  of  the  real  speech  program  "PROCESS."  N is 
the  number  of  processing  techniques  implemented  by  the 
system  (presently  twelve) ; M is  the  number  of  subprocess 
processing  steps  implemented  by  the  system  (presently 
twenty-four) . 

Figure  3 - Time  waveform  of  the  synthetic  vowel  /a/  as  generated  by 
"SPEECH";  (a)  no  noise,  (b)  signal-to-noise  ratio  of  0 db. 

Figure  4 - The  INTEL  process  applied  to  the  synthetically  generated 

vowel  /a/.  Each  waveform  shown  is  after  the  corresponding 
step  listed  in  the  first  paragraph  of  Section  V-A;  left  is 
for  a S/N  = «>  , right  is  for  a S/N  = 0 dB. 

Figure  5 - The  spectral  subtraction  process  applied  to  the  synthetically 
generated  vowel  /a/.  Steps  a-c  are  the  same  as  Figure  4. 
Steps  d and  e are  illustrated  after  the  corresponding  step 
listed  in  the  second  paragraph  of  Section  V-B ; left  is 
for  a S/N  = <=°,  right  is  for  a S/N  = 0 dB. 

Figure  6 - The  minimum  mean  square  error  filtering  process  applied  to 
the  synthetically  generated  vowel  /a/.  Steps  a-c  are  the 
same  as  Figure  4.  Steps  d and  e are  illustrated  after  the 
corresponding  step  listed  in  the  second  paragraph  of 
Section  V-C;  left  is  for  a S/N  = « , right  is  for  a S/N  = 

0 dB. 

Figure  7 - The  process  based  upon  pitch  extraction  appli.  i-  tj  the 

synthetically  generated  vowel  /s/.  Steps  p.  . re  che  same 
as  Figure  4.  Steps  d and  e are  illustrated  a',  ter  <_he 
corresponding  step  listed  in  the  second  paragraph  of 
Section  V-D;  left  is  for  a S/N  = » , right  is  for  a S/N  = 

0 dB. 


I.  INTRODUCTION 


The  understanding  of  speech  contaminated  by  the  presence  of  noise  is 
an  important  consideration  in  many  practical  communication  situations. 
Consider  the  airplane  or  helicopter  pilot  in  a noisy  cockpit  attempting 
to  communicate  with  ground  based  personnel;  or  consider  the  worker  in 
the  environment  of  noisy  machinery  attempting  to  communicate  via  tele- 
phone or  other  means  with  personnel  outside  the  noise  environment. 

There  are  two  cases  for  the  speech  in  noise  situation.  First,  there 
is  the  case  where  the  noise  is  present  at  the  speaker  (1-6);  second, 
there  is  the  case  where  the  noise  is  present  at  the  listener  (7-19) . 

Each  case  presents  a somewhat  different  situation.  In  the  case  of  noise 
at  the  speaker,  there  is  an  opportunity  to  suppress  the  noise  (relative 
to  the  speech)  prior  to  its  reception  by  the  listener,  thus  resulting  in 
an  enhanced  signal-to-noise  ratio  and  hopefully  more  intelligible  speech. 
In  the  case  of  noise  at  the  listener,  while  there  is  no  opportunity  to 
suppress  the  noise,  there  is  an  opportunity  to  process  the  speech  prior 
to  its  encounter  with  the  noise.  Conversely,  in  the  first  case,  since 
the  noise  is  already  present  with  the  speech  there  is  no  opportunity  to 
process  the  uncorrupted  speech  signal.  In  the  second  case,  since  the 
noise  is  in  the  environment  of  the  listener,  there  is  no  opportunity  to 
suppress  the  noise  level. 

There  are  several  types  of  noise  which  may  contaminate  a speech 
signal.  These  include:  impulse  noise,  large  amplitude  sine-wave  or  sum 
of  sine-wave  noise,  a conflicting  speaker,  and  wideband  random  noise. 

The  enhancement  of  speech  intelligibility  in  the  presence  of  each  of 
these  noise  types  is  under  study  by  the  Rome  Air  Development  Center  (1- 
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3,  20).  It  is  the  purpose  of  the  present  study  to  consider  the  enhance- 
ment of  speech  intelligibility  in  the  piesence  of  wideband  random  noise. 
A.  SPEECH  IN  NOISE  AT  THE  LISTENER 

One  obvious  method  for  improving  the  intelligibility  of  speech 


in  noise  at  the  listener  is  to  simply  decrease  the  level  of  the  noise. 
Although  there  exist  methods  for  decreasing  the  noise  levels  produced  by 


noisy  equipment,  these  methods  often  do  not  reduce  the  noise  levels 
sufficiently  or  are  too  expensive  or  too  inconvenient  to  be  practical. 

Another  obvious  method  for  improving  the  intelligibility  of 
speech  in  noise  at  the  listener  is  to  simply  increase  the  power  of  the 
speech  signal,  thus  resulting  in  a more  favorable  signal-to-noise  ratio. 
Although  such  a technique  may  work  in  certain  low-level  noise  situations; 
in  high  noise  levels,  the  need  to  conform  within  a pre-established 
maximum  sound  level  may  prohibit  the  use  of  this  simple  technique.  As  a 
result,  a method  for  enhancing  the  intelligibility  of  speech  in  high 
noise  levels,  without  increasing  the  signal  power,  is  desirable. 

Several  studies  related  to  the  intelligibility  of  speech  in 


noise  (at  the  listener)  have  been  reported  (7-19).  In  these  studies, 
various  techniques  for  processing  speech  (prior  to  reception  by  the 
listener)  have  been  evaluated.  Three  particular  techniques  (8,  10,  11) 
have  been  recently  reported  and  have  been  shown  to  offer  intelligibility 
enhancement  for  speech  in  white  noise  at  various  signal-to-noise  ratios 
and  in  one  case  in  the  environment  of  recorded  power  generating  noise 
(9) . These  three  techniques  involve  processing  speech  by:  (a)  high-pass 
filtering  followed  by  infinite-amplitude  clipping  (8) , (b)  high-pass 
filtering  (11)  , and  (c)  high-pass  filtering  followed  by  rapid  amplitude 
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compression  (8-9) . At  a signal- to-noise  ratio  of  OdB  (noise  at  90  dB 
(re.  .0002  dy/cm2))  an  intelligibility  of  greater  than  twice  normal 
speech  is  obtained  by  processing  speech  by  any  of  these  three  methods. 

At  other  signal-to-noise  ratios  a similar  intelligibility  enhancement  is 
also  achieved. 

B.  SPEECH  IN  NOISE  AT  THE  SPEAKER 

One  technique  which  may  sometimes  be  successfully  applied  for 
enhancing  the  intelligibility  of  speech  in  the  environment  of  noise  at 
the  speaker  makes  use  of  a noise-canceling  or  close-speaking  microphone. 
In  some  cases  such  a microphone  may  not  be  sufficient  because  the  noise 
is  at  too  high  a level.  In  other  cases  it  may  not  be  possible  to  use 
such  a microphone  because  it’s  too  inconvenient  for  the  speaker,  because 
the  speaker  is  unable  to  wear  such  a microphone,  or  because  the  speech 
is  already  added  to  the  noise  when  it  is  available  for  listening.  As  a 
result,  a processing  technique  for  enhancing  the  intelligibility  of 
speech  contaminated  with  noise  is  desired.  Several  techniques  have  been 
investigated  for  the  case  of  noise  at  the  speaker.  These  techniques 
include  high-pass  filtering  followed  by  infinite  amplitude  clipping  (4) , 
a technique  known  as  "INTEL"  which  involves  gating  the  second-order 
spectrum  and  the  subsequent  retransformation  back  into  the  time-domain 
(1-2),  and  methods  based  upon  the  use  of  linear  prediction  (5-6).  While 
the  techniques  investigated  for  this  case  have  displayed  encouraging 
results,  a significant  enhancement  in  the  intelligibility  of  speech  in 
noise  at  the  speaker  has  not  been  achieved. 

C.  THIS  REPORT 

While  both  cases  for  speech  processing  in  noise  have  many 
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practical  applications  and  while  both  cases  are  of  considerable  interest, 
greater  application  and  interest  exists  within  the  Air  Force  for  the 
later  case  (i.e.,  speech  in  noise  at  the  speaker).  As  a result,  the 
remainder  of  this  report  and  the  work  performed  in  its  preparation  are 
directed  specifically  toward  this  later  case. 

The  research  work  documented  in  the  following  sections  of  this 
report  has  included  an  analysis  and  study  of  methods  for  processing 
speech  so  as  to  result  in  an  enhancement  of  its  intelligibility  in  the 
presence  of  noise  at  the  speaker.  In  this  direction  two  specific 
approaches  have  been  pursued.  The  first  approach  involves  the  develop- 
ment of  interactive  software  for  the  analysis  and  study  of  processing 
techniques  applied  to  synthetically  generated  voiced  speech.  For  this 
study,  interactive  programs  have  been  developed  for  the  Honeywell  6000 
accessed  via  a Tektronix  4002A  qrj  for  the  generation,  processing,  and 
graphic  display  of  synthetically  generated  voiced  speech  sounds.  This 
interactive  system  is  described  in  Section  III  of  this  report  and  the 
results  obtained  with  this  system  are  described  in  Section  V. 

The  second  approach  involves  the  development  of  a general 
purpose,  batch-oriented,  speech  processing  system  for  real  speech  and 
its  use  in  the  study  of  speech  in  noise  processing  techniques.  This 
system  involves  magnetic-tape  communication  between  the  DDP-116  data 
converter  and  the  Honeywell  6000  computer.  Programs  for  analog  to 
magnetic-tape  conversion,  previously  developed  by  Captain  Robert  Curtis, 
are  utilized.  Programs  for  the  input  and  output  of  this  magnetic-tape 
to/from  the  Honeywell  6000  computer  and  general  purpose  programs  for 
processing  the  data  contained  on  these  magnetic-tapes  are  described  in 
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Section  IV  of  this  report.  In  Section  V the  results  obtained  using  this 
system  for  processing  speech  contaminated  with  noise  are  described. 

II.  DISCUSSION 

The  most  successful  present  technique  for  enhancing  the  intelligi- 
bility of  speech  in  noise  at  the  listener  is  the  technique  known  as 
"INTEL"  (an  acronym  for  INTelligibility  Enhancement  by  Liftering) . This 
technique  has  been  developed  under  the  direction  of  the  RADC  (1-2) . 

While  only  a small  speech  intelligibility  gain  has  been  achieved  by  this 
method,  a significant  enhancement  of  the  listenability  of  speech  in 
noise  has  been  demonstrated. 

It  is  one  purpose  of  this  work  to  investigate  the  operation  of  the 
INTEL  technique  through  both  an  examination  of  the  results  of  processing 
synthetic  voiced  sounds  and  through  the  processing  of  real  speech.  The 
results  of  this  examination  are  presented  in  Section  V. 

One  problem  with  the  INTEL  technique  is  that  It  requires  four 
Fourier  Transformations,  two  forward  and  two  reverse.  Since  calculating 
a Fourier  Transform  is  computationally  time  consuming,  it  is  of  interest 
to  explore  methods  which  require  fewer  transforms.  As  a result,  three 
additional  approaches  to  speech  intelligibility  enhancement  are  explored 
in  this  work  with  preliminary  results  documented  in  later  sections  of 
this  report.  These  three  approaches  are:  (a)  spectral  subtraction,  (b) 
minimum  mean  square  error  (MSE)  filtering  (21)  , and  (c)  methods  based 
upon  pitch  tracking. 

Spectral  subtraction  is  of  interest  since  it  is  computationally  and 
conceptually  relatively  simple.  It  involves  the  subtraction  of  the 
estimated  noise  spectrum  from  the  transform  of  the  speech  plus  noise 
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signal.  After  subtraction,  the  resulting  spectrum  is  retransformed  into 
the  time-domain.  Support  for  this  technique  can  be  gained  from  an 
intuitive  analysis  of  speech  plus  noise  and  from  its  seeming  similarity 
to  "gating"  which  is  performed  on  the  second-order  spectrum  in  the  INTEL 
technique.  One  problem  with  subtraction  is  that  it  requires  some  determina- 
tion of  both  the  magnitude  and  the  spectrum  of  the  noise.  In  the  INTEL 
technique,  the  magnitude  of  the  noise  is  automatically  accounted  for  (when 
gating  to  zero  is  used)  and  the  spectrum  of  the  noise  is  assumed  flat. 

In  the  results  presented  in  Section  V for  spectral  subtraction  it  is 
assumed  (as  a first  approximation)  that  the  noise  has  a flat  frequency 
spectrum  (i.e.,  white)  with  a magnitude  equal  to  the  average  first-order 
spectrum  magnitude  above  2.5  KHz. 

A method  for  determining  a filter  which  minimizes  the  mean  square 
error  between  a signal  embedded  in  noise  and  its  estimated  value,  assuming 
the  signal  and  noise  are  uncorrelated,  can  be  shown  to  be  given  by 
(21): 

Sss(w) (eq.  1) 

H(jw)  = Sss(w)+Snn(w) 

where : 

Sss(w)  = the  estimated  spectrum  of  the  signal,  and 
Snn(w)  = the  estimated  spectrum  of  the  noise. 

Using  the  longtime  average  speech  spectrum  for  Sss(w)  and  the  measured 
spectrum  of  the  noise  for  Snn(w) , the  transfer  function  H(jw)  can  be 
calculated.  Then,  multiplying  H(jw)  by  the  spectrum  of  the  input 
speech  signal  and  performing  an  inverse  transformation  results  in  the 


processed  speech  output. 


This  method  is  attractive  for  two  reasons.  First,  it  is  an  optimum 
method  (in  the  least  mean  square  sense)  for  separating  a signal  from 
noise.  Second,  it  is  computationally  simple  to  implement  since  it  is  a 
simple  filtering  process  which  can  be  performed  in  the  frequency  domain 
after  transformation,  or  perhaps  on  the  time-domain  signal  using  digital 
filtering  techniques.  A complication  of  the  method  is  that  it  requires 
an  estimate  of  the  noise  spectrum.  It  is  unclear  at  this  point  how  the 
results  of  this  method  vary  as  a function  of  the  error  in  estimating  the 
noise  spectrum.  In  the  experiments  described  in  Section  V it  is  assumed 
that  the  noise  spectrum  is  flat  with  an  average  magnitude  equal  to  the 
magnitude  above  2.5  KHz. 

Speech  in  noise  enhancement  through  methods  based  upon  pitch 
extraction  seem  intuitively  attractive.  This  is  a result  of  the  fact 
that  the  energy  of  a speech  signal  exists  at  harmonics  of  the  pitch 
frequency  while  the  energy  of  noise  is  distributed  throughout  the 
spectrum.  From  a knowledge  of  the  pitch  frequency,  those  lines  in  the 
spectrum  not  at  harmonics  of  the  pitch  frequency  can  be  suppressed 
leaving,  hopefully,  speech  enhanced  in  the  presence  of  noise.  Several 
problems  influence  this  method.  There  is  the  problem  of  accurately 
tracking  the  pitch  frequency,  even  in  a non-noise  environment.  There 
is,  on  the  one  hand,  a desire  to  analyze  speech  over  a long  time  segment 
to  gain  as  much  information  about  the  signal  as  possible.  There  is,  on 
the  other  hand,  a desire  to  analyze  a short  interval  so  that  changes  in 
the  pitch  frequency  during  the  analysis  interval  will  not  be  signi- 
ficant. A recent  report  by  Parsons  and  Weiss  (3)  suggests  that  the 
optimum  segment  size  should  be  40.8  msec.  A segment  of  51.2  msec  is 
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used  in  this  work. 
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Previous  results  obtained  with  methods  based  upon  pitch  tracking 
have  not  been  particularly  encouraging  (1-2,  20).  It  has  been  suggested 
that  even  when  the  pitch  frequency  is  accurately  determined,  for  example 
from  the  original  uncorrupted  speech  signal,  such  methods  have  not  been 
shown  capable  of  significant  speech  enhancement  in  noise.  The  use  of 

j 

comb  filtering,  for  example,  has  not  been  found  effective  for  improving 
the  intelligibility  of  speech  in  noise.  A discussion  of  the  use  of  comb 
filtering  for  speech  in  noise  enhancement  is  contained  in  a recent 
report  by  Weiss,  Aschkenasy,  and  Parsons  (1).  Another  report,  by  Weiss 
and  Aschkenasy  (2) , discusses  an  experiment  with  pitch  tracking  in 
which  good  results  were  reported  when  the  pitch  frequency  was  found 
-adequate.  Details  of  this  method  are  somewhat  sketchy.  However,  pointed 
out  in  this  report  are  the  difficulties  in  determining  pitch,  particularly 
at  low  signal-to-noise  ration  and  the  distortion  produced  by  inaccurate 
pitch  tracking  and  by  analyzing  unvoiced  speech  by  a harmonic  analysis. 

An  interesting  method  for  pitch  tracking  and  some  interesting 
results  obtained,  from  a relatively  crude  processing  method  based  upon 
pitch  tracking,  are  described  in  Section  V. 

III.  COMPUTER  SYSTEM  FOR  SYNTHETIC  SPEECH  PROCESSING 
A.  INTRODUCTION 

In  order  to  provide  some  insight  into  the  speech  in  noise 
situation  and  into  processing  speech  in  noise,  a terminal  interactive 
speech  processing  system  was  developed.  This  system  makes  use  of  a 
Tektronix  4002A  CRT  character/graphics  computer  terminal  connected  via 
telephone  modem  with  the  Honeywell  6000  Computer  System  at  RADC.  A 
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program  developed  by  Mr.  David  Clark  of  RADC  which  provides  FFT,  IFFT, 
data  input/output,  and  graphic  capability  was  used  as  a starting  point. 
Added  to  this  program  is  the  capability  for  analyzing  speech  by  the  four 
methods  (with  several  variations  each)  discussed  in  the  last  section. 

In  addition,  a program  for  generating  synthetic  voiced  speech  in  a 
format  acceptable  to  the  analysis  program  was  developed.  These  two 
programs  are  described  in  the  next  two  subsections  with  program  listings 
in  Appendix  A and  Appendix  B. 

B.  SYNTHETIC  SPEECH  GENERATION  PROGRAM 

A program  for  the  generation  of  synthetic  voiced  speech  called 
"SPEECH,"  is  listed  in  Appendix  A.  This  program  assumes  a sum  of  three 
damped  sine-waves  model  for  voiced  speech.  It  requests  the  following 
input:  (a)  amplitude,  frequency  (Hz),  and  damping  rate  (Hz)  for  three 
formants;  (b)  pitch  period  (msec);  and  (c)  RMS  signal  to  RMS  noise 
ratio.  The  program  generates  51.2  msec  of  the  speech  signal  (non-pitch 
synchronous,  with  sample  values  spaced  at  100  us),  adds  random  noise  of 
the  specified  signal-to-noise  ratio,  and  provides  the  capability  for 
viewing  any  generated  sample  values.  Finally,  the  program  outputs  the 
generated  speech  plus  noise  samples  into  a (previously  defined)  random 
file,  with  name  specified  in  the  RUN  statement  (program  line  number  0010 
in  program  listing,  Appendix  A).  The  format  of  this  file  is  acceptable 
as  input  to  the  synthetic  speech  processing  program.  The  program  runs 
under  GC0S  FORT, 

A sample  run  of  this  program  is  as  follows: 
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★ RUN 

nPJT*  AVPL.FREO, ALPHA,  FOR  ALL 
=1 .,730. ,27. 

=.b, IOVO. ,28.5 
=.04, 2240. ,46. 5 
I IP.  T PITCH  PRRIOD  IN  MSEC 
=2. 

NO  ICE? ( I -YES,  0-NO) 

= 1 

I 'p~T  S( RMS )/N( RMS ) IN  OB 

= 0. 

V I Ei»  ANY  SAMPLES? ( I -YES,  0-NO) 

= 0 


3 FORMANTS 


C.  SYNTHETIC  SPEECH  ANALYSIS  PROGRAM. 

A very  flexible  program  for  the  analysis  and  display  of  general 
time-waveforms  has  been  developed  by  Mr.  David  Clark  of  RADC.  This 
program  provides  the  capability  for  the  input  and  output  of  data  files, 
for  FFT  and  IFFT  analysis,  for  the  display  of  data,  for  data  generation 
under  cursor  control,  and  for  several  mathematical  operations  between 
two  data  files.  This  program  has  been  modified  to  include  the  processing 
of  speech  by  several  methods  based  upon  (a)  INTEL,  (b)  spectral  subtraction, 
(c)  minimum  mean  square  error  filtering,  and  (d)  pitch  extraction.  In 
addition,  the  updated  program  permits  the  capability  for  examining 
(plotting)  the  results  of  any  intermediate  subprocess  within  each  processing 
algorithm.  The  program  runs  under  GC0S  TF0RT,  is  named  M2  (following  an 
appellation  Ml  for  the  previous  program) , and  contains  instructions  upon 
initiating  a run  of  the  program.  A listing  of  M2  is  contained  in 
Appendix  B. 


IV.  COMPUTER  SYSTEM  FOR  REAL  SPEECH  PROCESSING 


A.  INTRODUCTION 

While  the  use  of  synthetic  voiced  speech  is  very  useful  for 
providing  an  understanding  of  both  the  general  speech  in  noise  situation 
as  well  as  specific  speech  processing  methods;  the  real  proof  of  a speech 
in  noise  enhancement  technique  is  how  well  it  performs  for  real  speech. 

In  an  effort  to  study  speech  processing  methods  on  real  speech  a set  of 
programs  were  developed  for  the  Honeywell  6000.  These  programs  run 
under  the  GCOS  operating  system  to  facilitate  the  input  of  data,  the 
processing  of  that  data,  and  the  output  of  results.  Since  the  Honeywell 
6000  does  not  provide  A/D  or  D/A  capability,  the  DDP-116  data  converter 
is  used.  Using  the  DDP-116,  input  speech  waveforms  are  sampled  (A/D 
converter)  and  these  samples  are  written  onto  magnetic  tape.  Having  a 
magnetic  tape  of  samples  to  be  output  to  the  DDP-116,  a time  signal  can 
be  constructed  using  the  D/A  converter.  Programs  for  the  DDP-116  for 
A/D  conversion  to  magnetic  tape  and  for  magnetic  tape  to  D/A  conversion, 
written  by  Captain  Robert  Curtis  of  RADC,  were  used.  Making  use  of 
these  magnetic  tapes  of  speech  data,  programs  have  been  written  for  the 
Honeywell  6000  Computer  System  to  read-in  these  magnetic  tapes,  process 
the  data  they  contain,  and  output  the  results  on  a second  magnetic  tape 
in  a format  suitable  for  input  (and  subsequent  D/A  output)  by  the  DDP-116. 

The  software  system  developed  for  the  Honeywell  6000  has  been 
purposely  written  in  a very  flexible,  very  general,  and  well  commented 
way  so  as  to  facilitate  its  future  use  at  the  RADC  for  any  general 
speech  or  signal  processing  task. 


B.  SPEECH  INPUT/OUTPUT  WITH  THE  DDP-116 


Input  signals  to  the  DDP-116  are  sampled  at  a 10  KHz  rate  (10- 
bit  sign-magnitude)  and  written  to  magnetic  tape  in  records  of  1024 
samples  (102.4  ms/record) . The  DDP-116  generates  three  tape  words  per 
sample  in  a tape  format  established  by  the  hardware  of  the  system. 

For  output,  the  software  of  the  DDP-116  requires  magnetic  tapes 
of  1024  samples/record  written  in  the  proper  sign-magnitude  format.  The 
output  software  created  for  the  DDP-116  can  read  these  magnetic  tapes 
and  generate  a time-waveform  from  the  sample  values.  The  D/A  converter 
is  10-bit,  sign-magnitude. 

For  both  input  and  output  the  DDP-116  requires  a one  word  zero 
record  at  the  beginning  of  any  magnetic  tape. 

C.  DDP-116 /HONEYWELL  6000  COMMUNICATIONS 

The  media  for  communication  between  the  DDP-116  data  converter 
and  the  Honeywell  6000  computer  is  7-track  magnetic  tape.  The  Honeywell 
6000  has  a memory  word  size  of  36-bits,  utilizes  two's  complement 
arithmetic,  and  when  communicating  with  a magnetic  tape  unit  utilizes 
six  tape-words/CPU  word.  The  DDP-116  has  a 16-bit  word,  utilizes  sign/ 
magnitude  arithmetic, and  communicates  with  a magnetic  tape  unit  with 
three  tape-words/CPU  word.  A method  for  tape  to  data-word  and  data-word 
to  tape  translation  was  developed  for  the  Honeywell  6000.  The  essence 
of  the  processing  necessary  for  this  magnetic  tape  communication  situa- 
tion is  illustrated  in  Figure  1.  Each  successive  word  read  from  magnetic 
tape  by  the  Honeywell  6000  consists  of  two  successive  words  written  by 
the  DDP-116  (6  tape-words/CPU  word  for  the  Honeywell  6000  and  3 tape- 
words/CPU  word  for  the  DDP-116) . Figure  1 illustrates  (right)  the 
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translation  which  must  be  performed  by  the  Honeywell  6000  to  map  a tape- 
read  word  into  two  data-words  and  vice-versa.  This  translation  process 
requires  a sequence  of  two  events.  First,  tape-word  bits  1-11,  14-17, 
19-29,  and  32-35  must  be  mapped  into  the  proper  bit  positions  in  the  two 
data-words  as  illustrated  in  Figure  1.  Second,  bit  0 and  bit  18  of  the 
tape-word  must  be  tested  to  ascertain  the  correct  polarity  of  the  data- 
word.  If  the  data-word  is  to  be  positive,  no  action  is  taken;  if  the 
data-word  is  to  be  negative,  the  corresponding  data-word  is  negated 
(two's  complement  sense). 

To  translate  a pair  of  data-words  into  a tape-word  requires  the 
reverse  of  this  process. 

The  subroutines  which  are  used  to  read  and  write  a tape  record 
by  the  Honeywell  6000  have  been  written  in  assembly  language  (GMAP) . 
These  programs  (called  RTB , WTB , and  WTBZ)  were  written  in  consultation 
with  Mr.  Albert  Proctor  of  RADC  and  are  included  in  Appendices  C,  D,  and 
E.  The  processes  used  to  translate  data-words  to  tape-words  and  tape- 
words  to  data-words  are  written  as  FORTRAN  subroutines  (called  DATTAP 
and  TAPDAT  respectively)  and  are  part  of  the  program  called  PROCESS 
described  in  the  next  section  and  listed  in  Appendix  H. 

One  additional  problem  exists  with  magnetic  tape  communication 
between  the  DDP-116  and  the  Honeywell  6000.  The  DDP-116  outputs  a 
parity  word  at  the  end  of  each  record  which  is  inconsistent  with  that 
expected  by  the  Honeywell  6000.  As  a result,  when  reading  DDP-116 
generated  magnetic  tapes  on  the  Honeywell  6000,  a parity  error  signal  to 
the  operator's  console  results.  This  signal  requires  acknowlegement  by 
the  computer  operator  in  order  for  the  program  to  proceed.  When  reading 
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a large  number  of  records  by  this  process,  the  operator  tends  to  prefer 
to  abort  the  job  rather  than  acknowledge  every  parity  error  (as  many  as 
2,000  during  a long  run). 

There  exists  a software  system  on  the  Honeywell  6000,  called 
UTILITY,  which  may  be  used  to  read  magnetic  tapes  ignoring  parity 
errors.  Using  this  system  a magnetic  tape  can  be  copied  onto  another 
magnetic  tape  (with  proper  parity  for  reading)  or  can  be  written  directly 
into  a data  file.  In  the  present  work  a copy  to  a second  magnetic  tape 
is  utilized.  This  second  magnetic  tape  can  then  be  left  as  part  of  the 
Honeywell  6000  system  and  the  program  "UTILITY"  need  not  be  used  again. 

The  program  UTILITY  runs  under  GCOS  CARDIN.  A listing  of  a 
batchjob  for  a tape-to-tape  copy  using  utility  is  shown  in  Appendix  F. 

A second  batch  program  to  read  a "UTILITY"  generated  magnetic  tape  into 
a data  file  is  shown  in  Appendix  G. 

D.  HONEYWELL  6000  SPEECH  PROCESSING  PROGRAM 

The  main  program  for  processing  real  speech  on  the  Honeywell 
6000  computer  system  is  called  PROCESS  and  is  listed  in  Appendix  H. 

This  program  runs  under  GCOS  CARDIN.  It  reads-in  speech  data  from  a 
file  (defined  in  the  batchjob),  processes  that  data,  and  outputs  the 
results  to  a magnetic  tape  (also  defined  in  the  batchjob)  in  a format 
acceptable  for  input  by  the  DDP-116.  The  processing  program  assumes 
that  some  method  of  overlapping  time-windows  (of  1024  data  samples) 
is  to  be  applied  to  the  input  data.  In  the  present  implementation  of 
the  program  only  a triangular  time-window  is  implemented.  This  processing 
program  is  purposely  written  very  flexibly  to  permit  its  future  use  for 
speech  in  noise  work,  for  other  speech  processing  work,  and  for  other 
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signal  processing  work. 

The  program  has  four  hierarchal  levels  of  subroutine  program 
depth  as  illustrated  in  Figure  2.  At  the  highest  level,  the  main 
program  performs  data  input,  data  output,  overlapping  time  windowing, 
output  amplitude  normalization,  and  calls  to  the  various  processing 
techniques.  At  the  second  level,  the  program  includes  the  processing 
techniques.  At  the  third  level,  the  program  includes  the  subprocessing 
steps.  The  fourth  level  (not  illustrated  in  Figure  2)  includes  sub- 
routines called  by  the  subprocessing  steps.  The  particular  processing 
techniques  to  be  called  may  be  inserted  with  appropriate  subroutine 
calls  at  program  lines  310  to  436.  A program  run  may  include  the  performance 
of  any  number  of  processes  (defined  by  NPROC,  line  102,  maximum  value  of 
20),  any  number  of  records  to  be  processed  by  each  process  (defined  by 
NREC,  line  105),  starting  with  any  record  (defined  by  NFREC,  line  107). 

The  program  creates  a magnetic  tape  of  the  processed  output  with  each 
successive  processed  speech  output  separated  by  ten  zero  records  which 
result  in  a 1 sec  silence  gap  for  use  when  listening  or  recording.  In 
addition,  the  program  prints  a listing  showing  each  processing  technique 
executed,  the  numbers  of  the  first  and  last  tape  records  created,  and 
the  CPU  execution  time  required  for  each  processing  method. 

In  its  present  version,  the  program  includes  an  implementation 
of  eleven  processing  techniques.  It  also  includes  twenty-three  third 
level  subroutines  for  call  by  each  processing  technique.  The  program  is 
well  commented  to  permit  ease  in  modification  and  future  usage. 

V.  RESULTS 

An  analysis  of  several  speech  in  noise  enhancement  methods  has  been 
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FIGURE  2.  An  overview  of  the  real  speech  processing  program  "PROCESS."  “ is  the  number  of  processing 
techniques  implemented  by  the  system  (presently  twelve) , M is  the  number  of  subprocessing 
steps  implemented  by  the  system  (presently  twenty-four). 
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carried  out:  (1)  using  synthetically  generated  speech  and  the  programs 
discussed  in  Section  III  and  (2)  using  real  speech  and  the  procedures 
and  programs  described  in  Section  IV.  This  analysis  has  provided  a 
basic  insight  into  several  speech  in  noise  analysis  methods  and  into  the 
speech  in  noise  situation  in  general.  The  four  subsections  of  this 
section  describe  each  of  the  four  analysis  methods  introduced  in  Section 

11  and  contain  representative  output  from  the  synthetic  speech  generation 
(SPEECH)  and  synthetic  speech  analysis  (M2)  programs.  In  addition,  the 
results  obtained  with  real  speech  are  discussed. 

For  the  synthetic  speech  illustrations,  results  for  the  vowel 
/a/  have  been  selected  for  all  four  methods.  Formant  data  is  from 
Peterson  and  Barney  (22)*.  A pitch  period  of  9 ms  was  selected  and 
signal-to-noise  ratios  of  co  (no  noise)  and  OdB  were  used.  Figure  3 
illustrates  the  output  from  SPEECH  for  these  two  conditions. 

For  the  intelligibility  tests  with  real  speech  two  utterances 
are  used.  Both  were  spoken  by  Captain  Robert  A.  Curtis  and  each  is  about 

12  sec  in  duration.  The  first  utterance  contains  speech  as  follows: 

"Testing . . .one. . . two . . . three. . .four. . . five ...  six. . .seven. . . eight . . .nine. . . 

ten...  we  were  away  a year  ago."  The  second  utterance  contains,  "We  were 

away  a year  ago .. .testing. . .one. .. two .. .three ...  four .. .five. . .six. . . 

seven. . .eight. . .nine. . .ten. . . may  we  all  learn  a yellow  lion  roar .. .Hawaii." 

The  first  utterance  is  in  a signal-to-noise  ratio  of  approximately  +6  dB; 

the  second  is  in  a signal-to-noise  ratio  of  approximately  -6  dB.  While 

this  corpus  of  speech  data  is  hardly  sufficient  to  quantitatively  evaluate 

*F1:  ampl,=  1,  freq.  = 730  Hz,  alpha  = 27  Hz;  F2:  ampl.=  .5,  freq,=  1090  Hz, 
alpha  = 28.5  Hz;  F3:  ampl.=  .04,  freq.=  2240  Hz,  alpha  = 46.5  Hz 
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speech  intelligibility,  it  is  sufficient  to  make  a preliminary  study  of 
the  processing  methods  and  to  make  several  qualitative  observations . 


A.  INTEL  PROCESSING  METHOD 
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The  INTEL  algorithm  is  described  in  two  recent  reports  (1-2). 

This  author's  interpretation  of  the  algorithm  from  these  reports  indicates 
that  the  analysis  of  a single  time-frame  of  512  points  (51.2  ms)  con- 
sists of  the  following  sequence  of  steps: 


(a)  Input  512  time  samples; 

(b)  Apply  a triangular  window  to  the  time  samples; 

(c)  Perform  a 512  point  FFT ; 

(d)  Set  to  zero  the  magnitude  spectrum  above  2.5  KHz; 

(e)  Square  root  the  magnitude  spectrum; 

(f)  Reverse  the  signs  of  all  odd  numbered  magnitude  harmonics; 

(g)  Perform  a 512  point  FFT  on  the  magnitude  spectrum  as  a real 
signal  with  zero  imaginary  part; 

(h)  Set  the  magnitude  of  the  five  low  frequency  harmonics  to 
zero  (gating) ; 

(i)  Perform  a 512  point  IFFT; 

(j)  Reverse  the  signs  of  the  real  part  of  all  odd  numbered 
harmonics; 

(k)  Square  the  real  part  of  the  spectrum,  make  it  a magnitude, 
and  restore  the  phase  from  the  original  time  waveform;  and 

(l)  Perform  a 512  point  IFFT  to  result  in  the  output  time 


signal . 

The  input  time  waveform  and  the  result  after  each  processing  step  listed 


in  the  previous  paragraph  are  illustrated  in  Figure  4 for:  (1)  no  noise 
(left)  and  (2)  a signal- to-noise  ratio  of  OdB  (right) . 
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1.  Discussion 


A discussion  of  the  overhll  INTEL  speech  processing  method  can 
be  found  in  two  recent  reports  (1-2) . The  discussion  here  assumes  general 
familiarity  with  the  INTEL  method. 

From  Figures  4a  - 4c  the  effect  of  noise  on  both  the  time 
waveform  as  well  as  on  the  magnitude  and  phase  of  the  frequency  spectrum 
can  be  observed.  From  Figure  4c,  it  is  noticed  that  the  effect  of  noise 
on  the  frequency  spectrum  is  to  add  a randomized  baseline  shift  to  the 
magnitude  of  the  spectrum  and  a seeming  randomization  of  the  phase 
characteristic.  The  phase  of  the  spectrum  of  the  speech  plus  noise  signal 
is  apparently  not  without  information,  however,  since  greater  intelligibility 
has  been  achieved  with  INTEL  by  restoring  the  input  phase  rather  than 
discarding  it  in  producing  the  processed  output  (20). 

In  step  e of  the  INTEL  process  a square  root  of  the  magnitude 
spectrum  is  taken.  In  a recent  report  (2),  Weiss  and  Aschkenesy  go  to 
great  lengths  to  demonstrate  that  a justification  for  this  square  root 
operation  centers  around  the  fact  that  such  an  operation  causes  an 
increased  fraction  of  the  noise  energy  (relative  to  the  speech  energy) 
to  be  concentrated  in  the  region  near  the  origin  in  the  second  order 
spectrum  (Figure  4g)  than  would  otherwise  be  present.  It  is  therefore 
argued  (2)  that  such  a square  root  operation  makes  it  possible  to  more 
completely  remove  the  noise.  What  is  unclear  to  the  present  author  is 
whether  the  "extra  noise"  energy  found  near  the  origin  in  the  second 
order  spectrum  is  due  to  the  original  noise  being  moved  or,  more  simply, 
the  first  order  spectrum  flattening  which  is  caused  by  the  square  root 
operation.  Following  this  latter  argument,  it  is  then  not  a desirable 
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feature  of  square  rooting  which  causes  increased  energy  near  the  origin 
of  the  second  order  spectrum,  but  rather  a consequence  of  square  rooting. 

Steps  f and  j of  the  INTEL  process  both  involve  reversing  the 
signs  of  all  odd  numbered  harmonics.  The  only  effect  of  this  pair  of 
sign  reversals  is  to  reverse  the  scale  of  the  axis  of  the  second-order 

spectrum.  A computational  advantage  has  been  suggested  as  a justification 
in  a recent  report  (2). 

The  essence  of  the  INTEL  technique  is  observed  in  Figure  4g 
which  shows  the  second-order  spectrum  for  the  synthetic  vowel  /a/  in 
signal-to-noise  ratios  of  » and  OdB.  Note  the  large  peak  at  low 
"frequency"*  in  the  second  order  spectrum  for  the  signal  plus  noise.  As 
shown  in  Figure  4h,  this  peak  is  removed  by  "gating"  in  the  INTEL 
process. 

Figure  4i  illustrates  the  first-order  spectrum  after  "gating" 
and  Figure  4j  after  the  sign  reversals  are  "reversed".  Comparing  Figure 
4j  with  Figure  4e  (the  spectrum  of  the  signal  before  gating),  the  effect 
of  the  gating  operation  upon  the  spectrum  is  observed.  In  the  absence 
of  noise,  the  effect  appears  to  be  primarily  a shift  in  the  amplitude  of 
all  magnitude  harmonics.  In  the  presence  of  noise,  the  effect  appears 


*The  use  of  the  term  "frequency"  in  describing  the  second-order  spectrum  may 
be  confusing.  While  it  is  true  that  the  units  of  the  horizontal  axis  for 
the  second-order  spectrum  are  not  Hz  when  carried  from  the  original  time 
waveform;  the  units  are  Hz  if  we  consider  the  first-order  spectral  signal 
like  a time  waveform  as  seems  to  be  a convenient  way  for  viewing  the  INTEL 
Technique.  In  any  case,  a clear  distinction  of  which  spectrum  is  being 
discussed  when  using  the  term  "frequency"  will  be  made. 


to  be  a removal  of  the  baseline  (Figure  4e)  by  a downward  shift  plus  an 
added  distortion  in  the  frequency  range  above  2.5  KHz.  A more  careful 
observation  of  Figure  4j , however,  reveals  a decrease  in  spectral  energy 
at  about  1.5  KHz  and  a marked  increase  above  about  2 KHz.  This  distor- 
tion was  observed  in  a recent  report  (2)  and  can  be  theoretically  shown 
to  be  a consequence  of  the  gating  operation  in  the  second  order  spectrum. 
A method  for  automatically  compensating  for  this  distortion  is  described 
in  a recent  report  (2).  Time  did  not  permit  the  addition  of  this 
compensation  algorithm  to  the  INTEL  procedure  programmed  in  the  present 
work. 

In  step  k of  the  INTEL  procedure,  the  spectrum  of  the  signal  is 
squared  prior  to  retransformation  back  into  the  time-domain.  The 
justification  for  this  step  is  an  attempt  to  compensate  for  the  square- 
root  operation  in  step  e(2).  This  author  wonders  why  it  is  thought 
necessary  to  compensate  for  the  square-root  operation.  There  is  little 
reason  to  suspect  that  attempts  to  preserve  the  original  frequency 
spectrum  are  necessarily  beneficial,  or  even  desirable.  The  now  classic 
work  of  Licklider  and  Pollack  (27),  as  well  as  work  of  others,  clearly 
shows  that  speech  can  withstand  severe  frequency  as  well  as  amplitude 
distortion  without  a significant  loss  of  intelligibility. 

To  test  the  need  for  the  square  operation,  a comparison  of 
speech  processed  by  the  INTEL  procedure,  both  with  and  without  the 
square  operation  were  examined  by  listening.  Speech  processed  without 
squaring  sounded  very  much  like  the  unprocessed  original.  This  indicates 
that  apparently  the  square  operation  is  essential  to  the  INTEL  process. 
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This  surprising  result  (surprising  to  this  author  anyway)  leaves  open 
several  unanswered  questions  regarding  the  use  of  the  square  root  and 
square  operations  in  the  INTEL  process. 

One  question  remaining  is  why  square  rooting  (and  subsequent  squaring) 
has  been  observed  to  improve  the  intelligibility  of  speech  in  noise  relative 
to  not  square  rooting.  One  argument  is  because  the  result  of  these  two 
operations  is  to  modify  the  output  frequency  spectrum  in  such  a way  as  to 
emphasize  lower  frequency  components  relative  to  higher  frequency  components. 
Since  lower  frequency  components  tend  to  best  survive  the  noise,  speech 
square  rooted  (followed  by  squaring)  would  be  expected  to  sound  better 
(enhanced  listenability)  in  noise.  Whether  this  process  results  in 
enhanced  intelligibility  remains  to  be  demonstrated. 

The  argument  of  the  previous  paragraph  would  suggest  that  the 
root  used  in  the  rooting  operation  (step  e)  should  be  dependent  upon  the 
signal-to-noise  ratio.  This  is,  in  fact,  what  has  been  observed  with 
the  INTEL  Process  (2,  20),  root  factors  of  one-third  and  one-half  being 
found  best  dependent  upon  the  signal-to-noise  ratio. 

One  final  observation  about  the  INTEL  technique  can  be  made  from 
Figure  4e.  From  this  figure  it  is  observed  that  the  INTEL  process  does 
not  preserve  the  triangular  time  weighting  applied  in  step  b.  This  may 
or  may  not  result  in  any  significant  distortion  depending  upon  how 
successive  triangularly  weighted  windows  sum.  This  particular  problem 
has  been  previously  observed  and  some  work  at  compensating  for  this 
distortion  has  been  pursued  (20) . 

2.  Tests  with  Real  Speech 
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From  listening  to  the  two  utterances  (described  at  the  beginning 
of  this  section)  processed  by  INTEL,  the  advantages  of  processing  speech 
in  noise  (by  INTEL)  are  apparent.  Clearly,  speech  processed  by  the 
INTEL  technique  "sounds  better."  This  has  been  described  as  enhanced 
"listenability."  While  there  is  some  question  as  to  whether  there  is  an 
enhanced  intelligibility  through  the  use  of  the  INTEL  technique,  there 
seems  no  question  as  to  a perceived  improvement  in  the  signal-to-noise 
ratio. 

A recent  communication  (28)  in  evaluating  the  INTEL  technique  in 
signal-to-noise  ratios  of  -5  dB,  0 dB,  +5  dB,  and  no-noise,  found  little 
intelligibility  gain  for  INTEL  processed  versus  unprocessed  speech. 

3.  An  INTEL  Extension 

Figure  4g,  the  second-order  spectrum  for  speech  in  no  noise  and 
at  a signal-to-noise  ratio  of  0 dB,  clearly  reveals  an  increased  amplitude 
near  the  origin  which  results  from  the  additive  noise.  In  the  INTEL 
technique  described  with  Figure  4 and  Section  V-A-l,  a suppression  of 
this  second-order  spectral  peak  to  zero  was  implemented  as  described  in 
two  recent  reports  (1,2).  Looking  at  Figure  4g  (no  noise),  however, 
reveals  that  even  for  normal  speech  this  peak  is  non-zero.  In  fact, 
from  an  observation  of  several  synthetically  generated  speech  sounds 
without  noise,  it  appears  that  the  peak  in  the  second-order  spectrum 
near  the  origin  is  usually  approximately  twice  the  amplitude  of  the  peak 
due  to  the  formant  frequency  envelope. 

Pursuing  this  observation,  an  experiment  was  carried  out  in 

\ 

which  the  INTEL  technique  (step  h)  was  changed  to  suppress  this  low 
frequency  peak  (and  other  low  frequency  harmonics)  by  a factor  which 
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causes  the  zero  frequency  peak  to  be  twice  the  amplitude  of  the  formant 
peak  in  the  processed  output.  This  is  accomplished  by  calculating  a 
factor  in  the  second-order  spectrum  given  by: 

X _ amplitude  of  the  zero  frequency  component 

amplitude  of  the  maximum  formant  envelope  component 

and  dividing  each  of  the  ten  low  frequency,  second-order  harmonics  by 

half  this  factor. 

Upon  listening  to  the  result  in  comparison  with  INTEL  a small 
improvement  in  intelligibility,  listenability , and  naturalness  seemed 
apparent. 

The  use  of  this  modification  to  the  INTEL  technique  has  advantages 
other  than  (perhaps)  enhanced  intelligibility.  First,  it  should  result 
in  less  distortion  of  the  type  shown  in  Figure  4 j (and  described  earlier) 
since  the  "gate"  is  less  severe.  (This  fact  was  not  verified  with 
synthetic  speech.)  Second,  it  results  in  an  overall  system  automatically 
compensated  for  a changing  signal-to-noise  ratio. 

A method  of  suppressing  the  zero  frequency  peak  in  the  second- 
order  spectrum  by  fixed  factors  has  been  previously  investigated  with 
INTEL  (20) . Results  have  indicated  that  a suppression  of  the  peak  by 
about  one-third  worked  best  at  a signal-to-noise  ratio  of  0 dB . From 
Figure  4g,  it  can  be  noted  that  this  is  in  very  close  agreement  with 
causing  the  zero-frequency  peak  to  be  twice  the  formant  peak  as  described 
in  the  previous  paragraph. 

B.  SPECTRAL  SUBTRACTION  METHOD 

As  previously  discussed,  one  problem  with  the  INTEL  technique  is 
that  it  requires  four  Fourier  Transformations.  A process  which  is 
intuitively  appealing,  appears  similar  in  function  to  "gating"  in  the 
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INTEL  technique,  and  requires  only  two  Fourier  Transformations  is  simple 
spectral  subtraction. 

The  basic  method  implemented  for  spectral  subtraction  consists 
of  the  following  sequence  of  five  steps: 

(a)  Input  512  time  samples; 

(b)  Apply  a triangular  window  to  the  time  samples; 

(c)  Perform  a 512  point  FFT; 

(d)  Estimate  the  average  noise  level  from  the  magnitude  of  the 
spectrum  above  2.5  KHz,  subtract  this  level  from  the  magnitude  spectrum 
and  zero  all  frequency  components  above  3 KHz;  and 

(e)  Perform  a 512  point  IFFT  to  result  in  the  output  time 

signal. 

The  input  time  waveform  and  the  result  from  steps  b and  c above 
are  the  same  as  for  the  INTEL  technique  and  are  illustrated  in  Figure  4. 
The  results  for  steps  d and  e are  illustrated  in  Figure  5 for:  (1)  no 
noise  (left)  and  (2)  a signal-to-noise  ratio  of  OdB  (right) . 

1.  Discussion 

One  problem  which  results  when  subtracting  from  the  magnitude 
spectrum  involves  the  action  to  be  taken  when  a difference  results  in  a 
magnitude  of  less  than  zero.  The  method  employed  in  this  work  for  the 
subtraction  of  a level  from  the  magnitude  of  the  spectrum  (step  d previous 
paragraph)  is  slightly  more  complex  than  simple  subtraction.  If  the 
value  of  the  magnitude  of  a component  is  greater  than  the  noise  level, 
the  noise  level  is  subtracted  from  the  magnitude  value.  If  the  value  of 
the  magnitude  of  a component  is  less  than  or  equal  to  the  noise  level, 
that  magnitude  is  divided  by  two. 


From  Figure  5d,  it  appears  that  the  simple  method  employed  to 
estimate  and  subtract  the  noise  is  not  totally  effective.  Comparing 
Figure  5d  with  Figure  4k,  the  output  spectrum  from  INTEL,  reveals  that 
INTEL  appears  to  do  a significantly  better  job  at  reducing  noise  than 
the  subtractive  method.  Examining  the  time  waveforms  for  the  output  from 
the  spectral  subtraction  method  (Figure  5e)  and  INTEL  (Figure  4e)  reveals 
that  the  spectral  subtraction  technique  does  not  modify  the  time  waveform 
as  severely  in  the  absence  of  noise.  Probably  a greater  amount  of 
subtraction  (greater  noise  suppression)  would  be  beneficial. 

2.  Tests  With  Real  Speech 

The  results  obtained  for  the  processing  of  the  real  speech 
utterances  described  in  the  first  paragraph  of  Section  V confirm  the 
expectations  of  the  previous  paragraph.  The  signal-to-noise  ratio  does 
not  sound  particularly  enhanced  over  that  obtained  by  simply  zeroing  the 
frequency  range  above  3 KHz  and  the  intelligibility  does  not  appear 
improved . 

As  a second  experiment  with  subtractive  noise  cancellation,  an 
additional  processing  step  was  added  between  steps  d and  e of  the 
subtractive  method.  This  additional  process  involves  the  use  of  a 
filter  to  emphasize  the  second  formant  frequency  range.  The  characteristic 
of  this  filter  was  chosen  to  have  a rising  slope  of  18  dB/octave  below  1 
KHz,  a passband  from  1 KHz  to  2 KHz,  and  a falling  slope  of  12  dB/octave 
above  2 KHz.  At  a signal-to-noise  ratio  of  -6  dB,  the  speech  with  high 
frequency  emphasis  sounds  less  intelligible  than  that  with  subtraction 
alone.  This  is  probably  due  to  the  fact  that  enhancing  the  second 
formant  frequency  range  is  detrimental  since  (at  this  low  signal-to- 


noise  ratio)  the  second  formant  range  is  so  heavily  obscured  by  noise. 

At  a signal-to-noise  ratio  of  +6  cB  it  is  not  clear  whether  the  high 
frequency  emphasis  is  helpful  or  not.  Certainly,  such  an  emphasis 
causes  the  speech  to  be  less  natural  sounding  in  both  cases. 

C.  MINIMUM  MEAN-SQUARE-ERROR  FILTERING 

As  indicated  in  Section  II  of  this  report,  the  method  of  minimum 
mean  square  error  filtering  is  an  attractive  technique  for  processing 
speech  in  noise  for  two  main  reasons.  First,  it  is  an  optimum  method  (in 
the  least  mean  square  error  sense)  for  filtering  a signal  in  noise;  and 
second,  it  can  be  implemented  in  a computationally  efficient  manner 
(relative  to  other  techniques) . 

The  implementation  used  for  this  method  is  based  upon  an  analysis 
from  Papoulis  (21)  as  described  in  Section  II  of  this  report.  The  implemen- 
tation consists  of  the  following  sequence  of  five  steps: 

(a) .  Input  512  time  samples; 

(b)  . Apply  a triangular  window  to  the  samples; 

(c)  . Perform  a 512  point  FFT; 

(d)  . Estimate  the  expected  noise  level  from  the  magnitude  of 
the  spectrum  above  2.5  KHz,  estimate  the  expected  signal  from  the  long- 
time average  for  normal  speech,  and  modify  the  magnitude  spectrum  using 
the  filter  H(jw)  as  given  by  Eq.  1 (Section  II);  and 

(e)  . Perform  a 512  point  IFFT  to  result  in  the  output  time 

signal. 

The  input  time  waveform  and  the  results  from  steps  b and  c 
above  are  the  same  as  those  illustrated  in  Figure  4.  The  results  for 
steps  d and  e are  illustrated  in  Figure  6 for:  (a)  no  noise  (left)  and 
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in  the  second  paragraph  of  Section  V-C;  left  is  for  a S/N 
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(2)  a signal-to-noise  ratio  of  0 dB  (right). 

1.  Discussion 

A discussion  of  the  minimum  mean  square  error  filtering 
process  is  contained  in  Section  II  of  this  report.  A theoretical 
justification  of  Eq.  1 can  be  found  in  Papoulis  (21) . 

From  Figure  6 it  can  be  observed  that  the  magnitude  of  the 
frequency  spectrum  of  the  ouput  (Figure  6d)  with  and  without  noise  are 
surprisingly  similar.  The  output  spectrum  after  noise  removal  (right) 
appears  very  free  from  noise.  The  output  time  waveform  after  noise 
removal  (Figure  6e)  also  appears  very  noise  free. 

One  additional  advantage  of  the  minimum  mean  square  error 
filtering  method  of  processing  speech  for  an  enhancement  in  noise  is 
that  the  method  is  directly  extendible  for  processing  speech  in  nonwhite 
noise  environments.  The  only  change  required  is  a modification  of 
Snn(w)  in  Eq.  1.  In  addition,  if  a technique  were  developed  for  detecting 
the  presence  of  speech  in  noise,  the  absence  of  speech  could  be  used  to 
continually  update  a running  noise  estimate.  This  running  noise  estimate 
could  then  be  readily  incorporated  into  Eq.  1. 

2.  Tests  With  Real  Speech 

Like  Figure  6 displays,  the  results  achieved  for  processing 
real  speech  in  noise  (signal-to-noise  ratios  of  +6  dB  and  -6  dB)  by  the 
method  of  minimum  mean  square  error  filtering  are  very  encouraging.  The 
signal-to-noise  ratio  is  significantly  enhanced,  the  naturalness  unchanged, 
and  the  intelligibility  sounds  improved. 

Two  experiments  were  carried  out  with  real  speech,  one  using 
the  procedure  outlined  in  the  second  paragraph  of  this  subsection  (and 
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illustrated  in  Figure  6) , the  other  with  the  addition  of  a zeroing  of 
all  magnitude  components  above  2.5  KHz  between  steps  d and  e.  In  both 
cases  the  results  were  very  good  with  little  observable  difference 
between  them.  An  examination  of  Figure  6d  indicates  that  there  is  very 
little  energy  above  2.5  KHz  such  that  little  change  would  be  expected. 

For  many  reasons,  the  method  of  minimum  mean  square  error 
filtering  appears  to  offer  a great  potential  for  speech  in  noise  intelligi- 
bility enhancement. 

D.  METHODS  BASED  UPON  PITCH  TRACKING 

As  described  in  Section  II  of  this  report,  methods  of  enhancing 
the  intelligibility  of  speech  in  noise  based  upon  pitch  analysis  are 
intuitively  attractive.  However,  as  also  indicated  in  Section  II,  such 
methods  have,  in  general,  produced  discouraging  results  (20). 

In  order  to  experiment  with  methods  based  upon  pitch  tracking, 
a technique  was  implemented  which  consists  of  the  following  sequence  of 
steps: 

(a)  Input  512  time  samples j 

(b)  Appl*y  a triangular  window  to  the  samples; 

(c)  Perform  a 512  point  FFT; 

(d)  Determine  the  pitch  frequency  using  a method  to  be  described 
and  zero  all  magnitude  components  between  pitch  harmonics;  and 

(e)  Perform  a 512  point  IFFT  to  result  in  the  output  time 

signal. 

The  input  time  waveform  and  the  results  from  steps  b and  c 
above  are  the  same  as  for  the  INTEL  technique  and  are  illustrated  in 
Figure  4.  The  results  after  steps  d and  e are  illustrated  in  Figure  7 
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for:  (1)  no  noise  (left)  and  (2)  a signal-to-noise  ratio  of 
0 dB  (right). 


1.  Discussion 

A discussion  of  the  use  of  pitch  extraction  methods  for 
speech  intelligibility  enhancement  in  noise  is  contained  in  Section  II 
of  this  report. 

The  method  used  for  estimating  the  pitch  frequency  from  the 
spectrum  of  speech  plus  noise  consists  of  determining  that  integer 
frequency  (FQ)  between  80  Hz  and  250  Hz  which  maximizes  the  following 
function: 


1 

N 


r 


i=l 


F(i-F0) 


(Eq.  2) 


where : 

N = the  greatest  integer  such  that  N'FQ  <_  3000  Hz 

F(i*Fc)  is  that  line  of  the  magnitude  spectrum  closest  to  the 
frequency  i’F0. 

The  results  of  Figure  7d  display  the  expected  spectrum  of 
the  output.  With  several  tests  on  synthetic  voiced  sounds,  at  several 
signal-to-noise  ratios  from  +6  dB  to  -6  dB,  the  pitch  frequency  was 
determined  quite  accurately  (to  within  a few  percent)  using  the  method. 

The  output  time  waveform  (Figure  7e)  indicates  the  appearance 
of  considerable  improvement  for  the  enhancement  of  speech  in  noise  by 
using  this  method  based  upon  pitch  tracking. 


M. 
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2.  Tests  With  Real  Speech 


Three  experiments  based  upon  pitch  extraction  and  the 
method  outlined  earlier  in  this  subsection  were  performed  with  real 
speech.  These  involve  three  degrees  of  suppression  of  non-pitch  magnitude 
components  in  step  d.  These  three  degrees  of  suppression  are:  (1) 
suppression  to  zero,  (2)  suppression  by  a factor  of  two,  and  (3)  suppression 
by  a factor  of  four. 

When  the  non-pitch  components  are  suppressed  to  zero,  a 
strong,  low  frequency  distortion,  probably  a result  of  the  time  analysis 
window,  is  apparent.  As  the  degree  of  suppression  is  lessoned  the 
strong  window  distortion  diminishes,  however,  at  the  expense  of  increased 
noise.  At  a high  suppression  of  non-pitch  harmonics  (suppression  to 
zero),  there  is  a noticable  loss  of  intelligibility,  particularly  at  a 
signal-to-noise  ratio  of  -6  dB . This  is  probably  due  to  the  inaccuracy 
in  pitch  tracking  which  is  greater  at  lower  signal-to-noise  ratios. 

With  a suppression  of  non-pitch  harmonics  to  other  than  zero,  the  method 
is  able  to  tolerate  greater  pitch  tracking  errors  without  as  serious  a 
degradation  of  speech  intelligibility. 

To  determine  some  measure  of  the  accuracy  of  the  pitch 
tracking  algorithm  applied  to  real  speech,  two  simple  experiments  were 
performed.  First,  a listing  was  created  of  the  pitch  values  determined 
for  the  utterance  at  a signal-to-noise  ratio  of  +6  dB.  While  no  standard 
of  comparison  for  the  determined  pitch  values  was  available,  the  printed 
values  seemed  reasonable  for  a male  speaker.  The  printed  values  were 
generally  fairly  continuous  with  fundamental  frequency  values  between 
110  Hz  and  120  Hz  during  voiced  speech  intervals. 


As  a second  experiment,  a constant  frequency  of  10  Hz  was 
subtracted  from  each  pitch  measurement  and  the  result  was  used  to  generate 
the  output  waveform  in  the  usual  manner.  Upon  listening  to  the  result 
it  was  found  unintelligible. 


Overall,  the  intelligibility  results  using  this  pitch 


tracking  method  sound  of  lesser  intelligibility  than  either  the  INTEL  or 
the  minimum  mean  square  error  filtering  methods.  There  are,  however,  a 
number  of  improvements  which  could  be  added  to  the  basic  method  and 
which  might  substantially  improve  this  first  attempt.  It  is  this 
author's  feeling  that  methods  based  upon  pitch  tracking  offer  the  potential 


to  result  in  an  enhancement  of  speech  intelligibility  in  noise  and  that 


such  methods  should  not  be  overlooked  because  of  past  discouraging  results. 


VI  CONCLUSIONS 

This  study  has  explored  several  methods  for  the  enhancement  of 
speech  intelligibility  in  the  presence  of  wideband  random  noise  at  the 
speaker.  This  exploration  has  involved  a study  of  the  effect  of  the 
methods  upon  synthetically  generated  voiced  sounds  in  noise  as  well  as 
on  real  speech  in  white  noise  at  two  signal-to-noise  ratios. 

Four  basic  methods  have  been  investigated:  (a)  INTEL,  (b)  spectral 
subtraction,  (c)  minimum  mean  square  error  filtering,  and  (d)  methods 
based  upon  pitch  tracking.  Several  variations  in  each  basic  method  have 
been  tested  and  numerous  experiments  with  speech  in  noise  have  been 
performed.  The  experiments  with  synthetic  speech  have  provided  a substantial 
insight  into  not  only  the  four  methods,  but  also  the  speech  in  noise 
situation  in  general.  Through  the  experiments  with  real  speech,  some 
qualitative  results  for  the  speech  processing  methods  have  been  presented. 
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It  appears  from  the  experiments  with  synthetic  and  real  speech  that 
all  four  methods  have  the  potential  to  result  in  an  enhancement  of  the 
intelligibility  of  speech  in  noise.  While  the  qualitative  intelligibility 
results  obtained  during  this  work  have  indicated  that,  of  the  methods 
tested,  the  INTEL  method  and  the  method  based  upon  minimum  mean  square  error 
filtering  seem  to  result  in  the  greatest  enhancement  of  speech  in  noise, 
the  other  methods  should  not  be  discarded.  There  are  several  reasons  to 
suspect  that  substantial  improvement  in  the  intelligibility  of  speech  in 
noise  can  be  obtained  through  the  use  of  each  of  the  four  basic 


techniques. 


'll.-,* 


REFERENCES 

1.  M.R.  Weiss,  E.  Aschkenasy,  and  T.  W.  Parsons,  "Study  and  development 

of  the  INTEL  Technique  for  Improving  Speech  Intelligibility,"  Technical 
Report  No.  RADC-TR-75-108,  Rome  Air  Development  Center,  Griffiss  Air 
Force  Base,  New  York,  April  1975.  (AD#  A011709) 

2.  M.  R.  Weiss  and  E.  Aschkenasy,  "Automatic  Detection  and  Enhancement  of 
Speech  Signals,"  Technical  Report  No.  RADC-TR-75-77 , Rome  Air  Develop- 
ment Center,  Griffiss  Air  Force  Base,  New  York,  March  1975.  (AD//  B004521L) 

3.  T.W.  Parsons  and  M.R.  Weiss,  "Enhancing/Intelligibility  of  Speech  in 
Noisy  or  Multi-talker  Environments,"  Technical  Report  No.  RADC-TR- 
75-155,  Rome  Air  Development  Center,  Griffiss  Air  Force  Base,  New 
York,  June  1975.  (AD#  A013767) 

4.  I.B.  Thomas  and  A.  Ravindran,  "Intelligibility  Enhancement  of  Already 
Noisey  Speech  Signals,"  J.  Aud.  Eng.  Soc.,  22,  234-236,  1974. 

5.  M.R.  Sambur  and  N.S.  Jayant,  "LPC  Analysis/Synthesis  from  Speech 
Inputs  Containing  Quantizing  Noise  or  Additive  White  Noise,"  IEEE 
Trans.  Acoustics,  Speech,  and  Signal  Processing,  ASSP-24,  488-494, 

1976. 

6.  S.F.  Boll,  "Improving  Linear  Prediction  Analysis  of  Noisy  Speech 
by  Predictive  Noise  Cancellation,"  1977  IEEE  International  Conf.  on 
Acoustics,  Speech,  and  Signal  Processing,  10-12,  1977. 

7.  R.J.  Niederjohn,  "A  Comparison  of  Three  Recently  Reported  Methods 
of  Processing  Speech  for  the  Enhancement  of  Speech  Intelligibility 
in  High  Noise  Levels,"  Proceedings  of  the  Midwest  Symposium  on 
Circuits  and  System,  325-331,  1976. 

8.  R.J.  Niederjohn  and  J.  H.  Grotelueschen,  "The  Ehancement  of  Speech 
Intelligibility  in  High  Noise  Levels  by  High-pass  Filtering 

Followed  by  Rapid  Amplitude  Compression,"  IEEE  Trans.  Acoust. , Speech, 
and  Sig.  Proc.,  ASSP-24,  277-282,  1976. 

9.  R.J.  Niederjohn  and  J.  H.  Grotelueschen,  "Speech  Intelligibility 
Enhancement  in  a Power  Generating  Noise  Environment,"  IEEE  Trans. 
Acoustics,  Speech,  and  Signal  Processing,  accepted  for  publication. 

10.  I.B.  Thomas  and  R.J.  Niederjohn,  "The  Intelligibility  of  Filtered 
Clipped  Speech  in  Noise,"  J.  Aud.  Eng.  Soc.,  18,  299-303,  1970. 

11.  I.B.  Thomas  and  W.J.  Ohley,  "intelligibility  Enhancement  Through 
Spectral  Weighting,"  1972  Conference  on  Speech  Communication  ani 
Processing,  360-363,  1972. 

12.  J.P.  Egan  and  F.M.  Wiener,  "On  the  Intelligibility  of  Bands  of 
Speech  in  Noise,"  J.  Acoust.  Soc.  Am.,  18,  435-441,  1946. 


43 


13.  I.  Pollack  and  J.M.  Pickett,  "Masking  of  Speech  by  Noise  at 
High  Sound  Levels,"  J.  Acoust.  Soc.  Am.,  30,  127-130,  1958. 

14.  I.  Pollack  and  J.M.  Pickett,  "Intelligibility  of  Peak-Clipped 
Speech  at  High  Noise  Levels,"  J.  Acoust.  Soc.,  Am.,  31,  14-16, 

1959. 

15.  J.C.R.  Licklider,  "Effects  of  Amplitude  Distortion  on  the 
Intelligibility  of  Speech,"  J.  Acoust.  Soc.,  Am.,  18,  429-434,  1946. 

16.  G.A.  Miller  and  S.  Mitchell,  "Effects  of  Distortion  on  the 
Intelligibility  of  Speech  at  High  Altitudes,"  J.  Acoust.  Soc. 

Am.,  19,  120-125,  1950. 

17.  W.  Wathen-Dunn  and  D.W.  Lipke,  "On  the  Power  Gained  by  Clipping 
Speech  in  the  Audio  Band,"  J.  Acoust.  Soc.  Am.,  30,  36-40,  1958 

18.  E.A.  Kretsinger  and  N.B.  Young,  "The  Use  of  Fast  Limiting  to 
Improve  the  Intelligibility  of  Speech  in  Noise,"  Speech 
Monogr.,  27,  63-69,  1960. 

19.  H.  Drucker,  "Speech  Processing  in  a High  Ambient  Noise  Environment," 
IEEE  Trans.  Audio  and  Electroacoust . , AU-16,  165-168,  1968. 

20.  Personal  communication  with  Dr.  Bruno  Beek  and  Captain  Robert 
Curtis,  Rome  Air  Development  Center,  Griff iss  Air  Force  Base, 

New  York. 

21.  A.Papoulis,  "Probability,  Random  Variables,  and  Stochastic  Processes," 
McGraw-Hill  Book  Company,  New  York,  1965,  pp  403-406. 

22.  G.E.  Peterson  and  H.L.  Barney,  "Control  Methods  Used  in  a Study 
of  Vowels,"  J.  Acoust.  Soc.  Am.  24,  175-184.  1952. 

23.  N.R.  French  and  J.C.  Steinberg,  "Factors  Governing  the  Intelligi- 
bility of  Speech  Sounds,"  J.  Acoust.  Soc.  Am.,  19,  90-119,  1947. 

24.  H.K.  Dunn  and  S.D.  White,  "Statistical  Measurements  on  Conversa- 
tional Speech,"  J.  Acoust.  Soc.  Am.,  11,  278-288,  1940. 

25.  I.B.  Thomas,  "The  Second  Formant  and  Speech  Intelligibility," 

Proc.  Nat.  Elect.  Conf.,  23,  544-548,  1967. 

26.  I.B.  Thomas,  "The  Influence  of  the  First  and  Second  Formants  on 
the  Intelligibility  of  Clipped  Speech,:  J.  Audio  Eng.  S.oc.,  16, 
182-185,  1968. 

27.  J.C.R.  Licklider  and  I.  Pollack,  "Effects  of  Differentiation, 
Integration,  and  Infinite  Peak  Clipping  Upon  the  Intelligibility 
of  Speech,"  J.  Acoust.  Soc.  Am.,  20,  42-51,  1948. 

28.  Personal  Communication  from  J.S.  Lim,  A.  V.  Oppenheim  and  L.D.  Braida. 

44 


•7 


< 

X 

M 

Q 

Z 

w 

04 

S! 


j 

m 


i 


u 

0) 

CO 

<u 

H 

60 

2 

< 

42 

y 

O 

(1) 

aj 

X 

CL 

LU 

co 

-a 

CO 

<u 

u 

CO 

-J 

*H 

4 

O 

-J 

^4. 

> 

< 

u 

♦ 

Or 

•H 

U 

co 

X 

Q) 

& 

4J 

c 

1ST 

. 

< 

>0 

X 

CO 

co 

X 

♦ 

^1 

co 

-< 

w 

w 

PL, 

C/5 


§ 

U 

60 

O 

M 

O. 

0) 


60 

c 


E 

n) 

U 

60 

o 

u 

o. 

c 

o 


cc 

< 

> 


M- 

r\) 

O 


cs 

LU 

cz 

u. 


a. 

< 


< - 


c 

< 


CJ 
L!  I 

»_i  ■ 

X 
C O 


— J 
* 


a 

it 


l-«  cn 
2 ro  * 
coa — 

•—  — 2.  II 

Y5  tJ-  *-i  *-i 

LU  . O 
roHM 
— • ii  rr 
Cmhc 
acta 
a. 


ro  — 

+ CO 

— r\i  * 

ii  • 

*-)  *— 4 

* w DC. 

x ■<  m- 

< > (M 

' > * O 

•— 4 CO  -K  — 

— »•— 1 a * 
a — x * — 
< n * *n 

> ■ • pj  <-> 

^ <\!  II 

• c ii  ^ o 

LU  C ^ •_! 

or  n r:  < C! 


M- 

LU 


\ 

H 

* 


tO 

* 

^r 

LU 


\ o 

H LU 
* tf> 

—> 

^ 2 
-J  *-• 


a.  -• 
x cz 
lu  a: 
* a. 

s 

— x 
*o 

2 H 

v->  *““* 

Dr-  X 
< 

cc  > H 

► + x 

— ^ x 
ii  •—  z 


C 

* 


LU 

4—4 

to 

4-4 

t-*  * 
LU  fT 

I— i fM 

X c 

co  — 


CM 

o 


^ + 
— 1 H 
— • IX. 


H U. 


X 

CO 


C IX. 
II 

C = C H 
O II  -Ii.  X 

M"  ' i— 4 • *— I 

4—44—1  2 CD  X 

*-  II  C CO 

< I—  CO  < Cr  X 

X CE 


— II 
II 


■ < X 
+ to  >-• 

*■">*  M | 


rn 

n 

2. 

4—4 

c 

to 

2 

y 

1 

X 

o 

2 

» 

S 

to 

4". 

LU 

CO 

>• 

2 

1 

Cr 

o — 

(M^  v 

— to 

— fM  O' 

CM 

Cl  — IU 

► H 1 

►UO  CO 

•XI  2 

— + — < 

an  x 

*-*  — • II  *-*  n c 


o -• 


o o 

<1  -o 


ii  cn 
o 


-•  c o < o < 


Cll  * X OJ  * CO 

- — ■ (—4  » . — . t — 4 » 

-4  2 C.  v 2 C, 

CO  4—4  *—4  <| 

Cl<  X-  LU  a.  cr  ixi 

x cc  — ■ x a 


x 

XI 

• II 

w ' 

II 

o 

X CM 
r-  x 
-c 

x c 

if)  cr: 


a co 

fM  co 

* 

o o o?coc  - c o 

— r\jcoM-'n'or^x>c 
COCOOOO'OC  — 

c c <*r  c e c c r-  c o 


o 


o 

xi 


■X 

•o 


o 

r- 


O 

fM 


■"!  c,  c c c c c o c c c o 

-fM^^X'OM)  > C — CM 

ccccoocooooc 


ccoocxccc 

n^  n'Ohx  >C  — 
r\i  ^ • *>.i  rvi  r\i  ''O  ^ 

CO  c c o o o o o 


A-l 


7~ 


z 

i 

o 


in 

iu 

>■ 

i 


f n 
< 


f\j 

x 

* 

m 

in 

co 

on 

an 

fM 

• 

o 

in 


o 

« * — 

CO 

<"X 

<M 

CM 

Or 

*~x 

u 

• 

\ 

w in 

*—* 

o 

fM 

r lu 

U, 

•— 

* 

fM 

— 

'O 

X _J 

• 

H 

* 

m 

C X 

•• 

p- 

<; 

* 

\ 

* 

2:  r-~ 

in 

/»■* 

*— 

«f“X 

H 

— 

o~s 

~j 

+ 

x < 

V.’ 

M 

UJ 

— 

in 

X 

Cl 

• 

* in 

Z 

V«* 

» 

• 

XI  >-• 

<*% 

v»* 

v— ' 

i*1 

c 

<* 

X 

fM 

» » 

X 

o 

< 

•< 

.— 

x > 

c 

X 

* 

Tl* 

< 

on  ^ — 

»w 

+ 

X 

\ rv j 

!T\ 

X 

X 

*— t 

* 

— 

< 

— “)  c 

-J 

-J 

10 

•—1  — 

< < 

* 

X 

-j 

• 

n m»  w 

CC 

< 

a. 

w 

x in  2- 

o 

— 

if 

* 

> 

»— 1 

N 

• 

*<  X 

o 

♦ 

i- 

H 

J5  * 

X 2 

.— 

X 

< 

V 

CK 

♦ 

'-<  O 

«*~N 

fM  « CJ 

«»~X 

O 

re 

-<  — 

+ LU 

•— 

• 

in 

II 

* 

X 

in 

II 

•— 

11  •n  < 

X 

• • 

• 

(X 

X II 

— , »— • 

* 

<3 

X 

-J 

>— 

n 

2’ 

< 

\ 

“>  > H 

rvj 

in 

m 

'.0  — < 

— ' > 

c 

C 

M* 

* 

o 

V— * 

CO 

< 

c 

_ 

— UJ 

1 

n 

n 

II 

v z 

V x 

w* 

H V O 

H 

X 

••> 

X— ^ 

o ~ c 

Q Q 

> 

—i 

-i 

_J  o 

< * 

• 

UJ 

< 

» 

c 

IU 

< 

/"■x 

rvj 

UJ 

fM  LU 

#— i 

22H 

ii 

Cl 

x on 

II  H 

o 

H 

— 

H 

r-» 

-J 

— 

H 

• 

-HJO. 

H < < x 

Q: 

y 

=?  fM 

*"s  Z 

< 

» — < 

CC 

< 

►— » 

CC 

_i 

m 

>-1  o 

H — j c 

n 

U x x o 

CC  23 

< 

< 

< 

t — < *— < 

LU 

'JL 

X 

UJ 

^■x 

Cr 

^•X 

^■x 

< 

ii 

D: 

II 

C x < H 

II  II 

in 

^X 

x»' 

H 

o 

<x 

Q. 

w or. 

X 

•—* 

e 

(JL  X Q 

.JT 

LL 

5 

o 

oq  Ci  r c 3 in  in 

— fM 

it 

LU  Z 

in 

'/) 

x 5 

< X 

u. 

X X >- 

X X UJ 

o 

O 

O 

r-» 

O 

rv 

on  n 

c 

c 

X 

•— 

fM 

fM 

fM  "0 

\ 

> 

— 

X 

*“ 

— 

e 

O 

c o 

o c 

C 

c 

o o o 

c 

O O 

r 

o 

o c o o c 

o c o o 

f-x 

fJN 

o 

fM 

00 

in 

■c  0“ 

00 

o 

o 

•— 

fM  on 

IflOMO 

> o 

— fM  on  m-  n 

<c  r-  to  x 

c 

fM 

nj 

on 

!~> 

-V*  *V| 

"O  ~i~ 

*VX 

r**s 

»+ 

t 

*-*- 

Tf 

«r^ 

*t 

•t 

r*- 

■n 

•n  n n on 

n 

n n n 

\r* 

o 

n 

1- 

O 

o 

O O 

o o 

o 

O 

c 

c c o 

o 

o o e c 

ooccoocoooo 

o c o 

u 

A -2 


iMnft  . HU  m*ti± 


Source  listing  of  *\«  program  "M2”,  synthetic  speech  processor. 
Only  additions  to  the  pre  :essor  of  this  program  "Ml"  (created  by  Mr. 


**•  ■ — 


O 

c 

in 


o 

c 


00 


H 

Cr. 

< 


•a 

<u 

■u 

to 


1) 

U 

« 


X. 

u 

(0 


TJ 
•H 
> 

o 


c 

H 

c 

X 

!M 

j. 

— 

X 

9 

r-* 

o 

-i 

o 

c 

Cl 

— 

_ 

’li 

» 

. 

C> 

m 

9 

LU 

* — 

o 

M 

. 

o 

z 

o 

(M 

— 

• 

► 

9 

9 

cr 

Cr 

oo 

X 

o 

■ — “ 

s . — 

o 

CO 

>- 

m to 

e — - 

o 

ip 

CO 

— 1— 

LU  00 

► 

UJ 

CD  LU 

x c — 

8 

o 

• 

o 

* -J 
~ CL 

z H ’ 

CO 

or 

oo  z 

z c — 

9 

CL 

fir 

► c 

o ► 

o 

CL 

oo  O 

H — O 

o 

OJ 

. — 

z o — 

r> 

9-9 

— CO 

X c 

t — « • * 

9 

— h— 1 

*r* 

1— « 

C C‘  cv 

o 

O H 

H 

X 

am  - 

c 

Dr 

^ — o c 

M 

9 Q- 

• 7) 

<5 

9 

(D 

C < H 

C cr 

A' 

Z Cu 
► L- 1 

CD  H ► 

Z Or  f- 

O 

KT. 

crinhC 

CO 

H X 

O 

UJ 

C c an  o 

J 

Q ■— 

H H < 

in 

o 

-(/ICC- 

< 

► OO  H 
> tn 

cun  9 

' i/A 

o 

o 

<3- 

9 

8 

00 

► 

o 

o 

OJ 

► 

o 

o 


H 

/“•* 

:o 

o 

•o 


c x o in  — 
cr  c ui  » uu 
any  ^ ji.o 
c x < — • 
j c a < a c 
ii'oarnr 
H —•ox 

— • iu  — > z a.  CL 
H cr  —>  in  x 
C Z Or  X 
H — • C C 
H *J  H .-1 
j J J _? 
i— )*-!_]<  _J  < 

< < _J  O -J  o 
O CJ  < < 

yooo 

O I C It. 
O I (M  I oo 

— I — I — 

- o o 

„ O C — O.' 

in  c ■<*  xi  in  x» 
X **  r\t  oj  -m  oo 
fMNOOOO 


I 

I 

» 

I 

O 


• • % 


O Ul  00  r\j  < 

Z H o — J 

— 1 z.  m 

(x.  — 1 >-  ^ cr 

x < < LU 
u i >jj  cr  H 
o?Q  nz 

C — ' < C LU 

at-  — • r 

CL  Z Z X •• 
^ ^ y f— < 
JCrSliiz 
LL.  23  » S-  — • 
<— i x C — 1 or 
z in  y c a 


H ~ ► 

x j. 

M -| 

cr  *>•  oo 

u>  H •• 
H cr  om 
!-  Z < -I 
x UJ  H — 

< s x ~ 

_i  « z, 

M (— < ••  O 

g z c. 

< — 1 < 

LU  Or  1X1  C 
Q-  CL  Or  CD 


I 

I 

O 

O O C n c 

O C <M 

n •c  c 'C  x 
o o c-  o o 


O — Oi  o 

00  oo  00  X/ 

■O  'T  'T  O 

o o o o 


B-l 


<"0 

D 

lu 

H 

< 

a 


a 


U.1 

X 

Ct: 

4 • 

CO 

• 

(X 

LU 

H 

'C 

CO 

in 

I.U 

G 

a 

o 

OJ 

Q 

H 

•».  • 

\ 

♦— • 

UI 

r-v 

UJ 

77. 

fM 

a 

7s 

— i 

Oj 

si' 

a 

a 

a 

>< 

< 

•S 

s — 

X 

i 

>- 

7 

>- 

S' 

n: 

• 

< 

OJ 

< 

a 

r- 

H 

a 

4—4 

m 

2T. 

»-4 

o 

a 

a 

H 

/■> 

CM 

LU 

(II 

in 

< 

CO 

— 

*w» 

CO 

CO 

a 

>- 

2: 

X 

w 

* 

< 

UI 

+ 

cr 

< 

C 

< 

c 

V 

• /•% 

"V! 

cr 

OJ 

ci 

Q. 

c 

a 

r> 

-o  — 

X 

CL 

■k 

UI 

r 

3 

a 

c? 

O 

in  •■ 

7* 

*-* 

77 

2 

* 

^ NJ 

tr.  < 

2r 

r~ 

^■N 

7 

zr 

7l  X 

T* 

fM  <■  O'  (M 

“ 

O' 

■— • 

a. 

7 

a 

a 

a 

' — 4 

Q: 

i 

a 

a 

“ Q: 

*— i 

\ in 

7} 

< 

H 

— 

o o 

rr  a 

X 

> 

77 

a. 

r> 

a 

m x co 

77 

*o  cm  i— 

jtj 

V/ 

| -4 

T- 

!-< 

f-  !- 

H1  c 

a 

< 

H 

< 

a 

*—• 

VHH 

X + 

uj  a z at  H 

UJ 

t > 

*— 1 

a 

ui  m 

a 

a 

4—4 

X 

OJ 

U) 

c 

X u:  7L 

LU 

* -)  C - 

*— 1 

o- 

C 

< 

w 

a a 

cr  a cc 

o. 

a 

o 

a 

» 

a 

X 

<— i c a 

r c 

r~s 

a.  a 

CM 

77 

oa  a 

~ < OJ 

^ jn 

< OJ 

a 

/>~s 

— _ 

a <->  3t 

•— ' \r\ 

— U. 

— 

(0 

— 

oj  tr. 

— 

— 

oj  ro  c 

— 

< 2J 

•— 

in  < > 

-c 

— 

r- 

•-  m 

X 

a 

t N LU 

H CM 

• II 

• 

77 

• 

!0 

• 

• 

U' 

in  H 

• cm 

• T O' 

• 

amt-* 

• 

‘O  ♦— « 

• 

V/  • 

» 

—>  *■" 

C 

< 

y 

o < 

• 

4—4 

4—4 

ry. 

Q:  O Ui 

• 

O oj 

c 

W' 

» 

cr 

r> 

o 

c 

Cl 

*/.  C X 

a:  — 

LU 

fv 

*— i LU 

W> 

UJ 

— 

c a 

• 

a a — 

• 

UJ 

— 

C U) 

► — * 

— u. 

LU 

2 

4—4  X 

< II 

'JL.  * 

• 

j— 4 

< 

• 

77 

II 

X 

X 

• 

X 

ii 

;n 

• /-k 

It 

/“N 

• a n 

o 

• 

H 

II 

CO 

• 

CO 

II  1 

• 

x 

< • 

a “•> 

II  C 

I—, 

—7 

H 

4— i 

♦— i 

II 

II 

[—4 

S-* 

»— i 

n 

H a 

♦— 1 

fi 

H a >-h 

ii 

H 

a. 

II 

H 

*-*  II 

H 

r_) 

X H 

m 

" * n 

maun 

CO 

X 

to  ~ 

co  a 

/"■s 

CO 

n 

CO 

a 

CO 

a 

O c r.  r. 

to  o 

— rj  < a 

<0.0 

— 

oj  < a O 

— 

< S’ 

O OJ 

< r>  c 

— 

< 

CO  o 

•— 

< 

to  o — 

< 

CO 

< c 

^ r— 

• + 

_i 

*— 4 

_1 

a 

CM 

a 

fO 

J >-H  •J- 

a _ tn 

» 

^1 

>o 

» 

a 

a 

f-  • 

a 

a a a 

< 

It 

CL 

X 

'W 

in 

4-^ 

♦— 1 

'w' 

a1 

♦— 1 

'W» 

♦—4 

*— 4 

'w' 

X 

*—4 

a 

►—4 

a 

> — ' ^ 

W s/ 

X 

< a > 

'w' 

X 

V 

X H 

3 

a a c 

s—* 

a 

a 

> 

c ^ 

11- 

v: 

< a H 

C"  X 

y 

X x 

4—4 

O 

•— i 

< n 

< < 

♦— 1 

< 

o 

X 

*->  iu  o a 

•-  a o a 

► — i 

< o a 

4—4 

:u  o a 

►-H 

< 

uj  *— 1 1 r 

H 

'U 

1 

co 

1 

Ls- 

I 

0^ 

i 

NJ 

1 

H 

| 

n 

1 

H 

1 

CO 

i 

1 

i 

> 

1 

1 

i 

i 

O 

1 

1 

c 

i 

i 

o 

1 

1 

O 

i 

C 

1 

o 

1 

1 

i 

i 

i — 

— 

1 

OJ 

i 

fO 

(M 

1 

i 

n 

1 c 

m 

1 

r-  o 

i 

x r~ 

1 

> i 

o 

cj 

CJ 

CJ 

CJ 

CJ 

CJ 

o 

CJ 

CJ 

c C 

r~ 

;*»  r*- 

c 

Q 

c 

o 

c 

c 

r~ 

N- 

o 

c 

c o c o c c.  c 

o 

c 

o 

c 

o 

c o c 

c 

C O C' 

•'t  m -o 

X) 

> 

C 

— 

oj  n 'j 

in  o N- 

x a c 

— OJ 

ro 

M- 

1^0  7 

X 

> 

c 

— 

OJ  oo 

't  m c 

»*s- 

X 

O'  C O' 

.o  -r 

\0  V* 

\T 

f- 

r- 

r-  r- 

r-  r- 

r- 

r-  no 

no  'O 

X' 

?0 

7 »Y>  ^ 

0* 

X 

O' 

"X 

O' 

\ 

CV  •> 

Os 

>> 

O'  n n 

occoccoocococoooo 

oocoocccc 

o o c c 

coco 

c 

o — — 

B-2 


< 

fM 

»»> 

•— 

• 

_ — 

fM 

» » 

7" 

zr 

— 

“ 

; 

O 

*— 1 *— 1 

cr 

X 

in 

U1 

CL  — 

— 

v-/  >-/ 

LL 

H 

n-x- 

H 

:u 

— 

< < 

LL 

*— 1 

0 

LL' 

A 0 

— 

'A 

• * 

•a 

Cr 

• 

^_j 

X H 

H 

*-s 

w- 

CL 

_ H 

z 

fM  fM 

< 

0 

fM 

X.  C Q: 

Cr 

UJ 

* * 

Cr 

X 

O < 

Q 

2F 

*— « *— • 

H 

< 

fM 

O 

H H 

ll 

LU 

— 

X O (A 

in 

-J 

< < 

Cx. 

0 

in 

f • 

>-*  • X. 

«< 

UJ 

w w 

Ul 

7 

' — 

'—4 

CO  * " 

«c. 

nj  <m 

►—4 

4-4 

LL 

a.  u j 4~  — 

Cl 

O 

•T  “» 

O 

LA 

Q- 

• rr  ^ 

H 

C 

< < 

X 

CA 

•z 

UJ 

OH  ^ 

‘ZZ. 

x c 

X 

*T" 

H H 

X.  0 ~ 

LU 

m 

1—4 

X Cr  r-  c 

act: 

cc 

cc 

CC 

cr 

< < cc  LL  - 

O X H 

(A 

— ■ < -0 

X UJ  0 

X 

X *J 

X 

X 

X c 

#*“ 

< O 

4""X 

HHC2 

H *— 1 — 

H — » 

H X 

H 

H 

~z  in  H uj  — 

Cr. 

y*' 

f\l  H 

cr  tA  ►*-1 

_ 

may 

UJ  (M 

LU  <•  CVJ 

UJ 

r\j  Uj 

►H  ^ 

1 ! » rS)  V 

CL 

*— « 

.(A 

< X LO  2 

>c 

o-  X x 

Or  — 

CC 

CC 

+ C: 

moaaz 

0 f\i  < 

H V-4  ' 

r\i  ^ C ^ r\j  — r\i  + * ^ LU  •-<  'JJ  ~ (A  ll  tj-  lu 

• c ll  — l — — e\jc  — — o — — — in  > ►_)  rn  x w m *—  - jc 

f)  — ~ lA  — ii">  * — '-■in  — — v-  lA  — — — X.  - >-<  O"  >-  CL  Or  •»  O',  •—« 

to  • u - • ll  • •;?  »’—'•>-  •*•111  ►**•>— iv^  o < < ui  u ■ [ — < » * — * 

r\j  o m 2 c - - o o o<  — >-  r>  n - —>— >o  x x lu  cc  i—  h cr  <m 

II  . • W Cr  *->  UJ  rr  II  • LL'  •—  II  LL  UJ  Or  II  ■»'  Ul  < ||  v*  LU  n *->  2 O'  2.  2 H r < ► OC 

o O • UJ  < • >—  — I o *KH  | • CX  >— • LL  • ~ LL  LL  • X.  < 2*-><CUJC/)UJt-*~< 

II  II  H > T H ll  H ll  H < II  H Cl  ll  ll  HCX  •—  H — < = < r (A  ~ _J 

^ tn  2 u in  c c « in  in  c ^ in  o ~ ir  o — ia  o o x n z m » j « 7.  x 

C - f\l  < H < CC  — <N  < (A  H — < LLJ  ^ — C:U'0(M  — < UJ  cc  xccxH~H~  CO 

_ » »_j  x _j  < — .._jcc  — .jot-  »ja-  * *_jia_jx  ottsuizo^.DHZ 

LI'-J'-’O  « LL!  •— ''-'C  «-■  HOLLKQ7r*-"<M<  < 

c^''U.y<iunc^u.>C''U.2C^u.Hr''<-Li.iJJ<iu2Li'3c*-'£rli|c:iiic*-< 
rvj.u.'-'^OM'-ciH'Lr.iiHoaiLHinciLi^.v'jcujQ.inoniiaaKOQ: 

He  rr  ir,  lu  < cr  H 

ol  I o i o I c T c V I I 

CO  I—  I r\j  — I rr  rn  I tJ-  *3-  IcO'C  I O I I 

— — I — I I I I — — I — I I 

CJ  O O O O OO  O 

C C re  0000c  coocccccocococcccccoccoooccec 
c r\jH  'tLT\'cr~x>c  — iMnnvin'Cr'Ti^in’cr®>c  - c — cm  m ^ in  c m o> 

*—  — — — — — — — — — — — — — — — — — — — — — — — — — — fMf\ifMfMfM<MfMfMfMfM 


2*30  F 0=  1 000 ./XI  'AX 

2*40  PRINT* "PITCn  FREQ.  ="  , F0,  " MSEC" 

2*50  IF ( LAST. EQ. 3 ) RETURN 

2 *60C CONSTRUCT  NEH  RE(F ) *IM(F)  KEEPING  ONLY  HARMONICS  AT  PITCH  FREQ. (6) 


2 170  6 DO  60  1 = 2.512 
218.1  F(  I,  1 )=0. 

24y'i  60  F ( I ,2 )=0 

2500  DO  70  I»IMA:<#3000.r«AX 


in  <-» 

• — CM  — 

+ * * n 

in  r--  — w 

cm  x x < 

— 7 7 II 

o x x ~ 

in  ^ — 

• — •<  •<  * 
> + II  II  x 

2'  <*■■>  o-s  ^0 

V =7  — r\)  Z. 
— ' 7'  * • X 
x J7r  I 

II  II  3 3^ 

>.  s r-  z — 

II 

**** 

X J JlL'l’L 


o c-  c c 

- (M  'f 

n »n  lOi.r'' 

<M  (V  CM  CM  CM  i 


II  C'ffw 
^ ^ l* J X 
<M  >C'  ’— 1 = 

. • cr 

-0  — 7’ 
33  Ul  C « i 

2.  • X < 
XH  X 
I i r.  u>  o 
’t  < (/) 


x < 

X 

U!  O S3 

in  or. 

X X 33 
X’  XH 
> < UJ 

2.  O Q: 


o o o < 

h cc  > 

^ in  o’ 

i cm  oj  cm  < 


O .l  i 

\ 

• cm 

CM 

s in  * 


H >-  II  < ■ 

CJ  X ■— i || 

x — < 
x in  o — c 
H ft  OC  * 

in  uj  >-i  ► 

c h n x •. 


<r> 

o 

coo 
ni  nr  Tt  i 
cr'C 

CM  CM  CM  i 


X 

*— 

CM 

* 

33 

S"" 

3-*. 

o 

X 

re 

► — * 

X 

•or 

2; 

3 

H* 

< 

u 

X 

X 

C 

y 

1 

i 

cn  ~ 

< 

in 

M-  2 

x.  = 

X 

• 

• — CM 

— 

— C 

<C  C7 

H 

o 

+ - - 

in  in  x 

O. 

70 

r~ 

in  y y 

H V 

3> 

C 

CM  33  33 

<* 

< UJ 

sr 

CO 

ro 

— Z 3~ 

II 

ii  a 

D*  ' 

• 

O XX 

/"■<* 

mu 

UJ 

X 

in  — — 

0— 

r\!  cr 

HH  3 

'0 

< 

. — < c 

• 

• • 

O'  W 

i — i 

5t‘ 

r>  + li  ii 

~ C X X 

«'*> 

• 

HH 

0—  JT  y~s 

rj  ui 

o • — * 

X. 

c 

II 

Nr-iv 

✓ • 

r".  • 

x < 

II 

l-H 

*-<  33  * » 

XXo 

X 

x 

x x y r 

1 

i in  uj  o 2 

< 

CM  O 

II  II  X 37 

m-  < in  x 

X 

• 

O *“H 

ss  r-  x.  r. 

— 

(TX3 

H 

l-H 

II 

x x x x n 

lO  ~ 

u:  X I-*  c. 

r \ 

V_/ 

73  " — N — ■> — X y 

I J IL  U.  1U1  M Z U 


C pc- 

3>  C. 

o -v. 

fM  CM  CM  CM 


O I 

a i 

o __ 

— i\  C0  if 

r—  r-  r-'- 
CM  OJ  CM  CM 


x rr  i 
a I'.'  x 
to 

I 

I 

I 

I 

o 

C O c 3 

m <o  c 

-»  r-  C 
CM  CM  O 


PI 

• 

o 

o J 

s 

•c 

LL 

• 

in  u.' 

Ct  — 

m 

»— « 

* 

*— • 

H 

UI 

rg 

H 

x-» 

o in 

33  C 

\ 

f\J 

O 

rg  x 

• _» 

•r  H 

^x 

.<• 

x CL 

33  H 

¥ 

a- 

• CM 

(M  > 

7 c a 

X 

H 

pg  O 

2 

* c~ 

O < 

i 

X-* 

r\» 

tr 

■K  ^ 

X 

C 

rg  O 

H x H 

• 

_i 

77 

* 

x-» 

H 

— 

37  o in 

r- 

Ui 

X 

C 

pg 

< 

m 

m C 

►H  • £=” 

in 

> 

W 

rg  H 

X-P 

H 

CO  - — 

rg 

111 

LL 

>- 

• 

-c 

X " 

IX. 

cl  uj  x _ 

* — 

J 

+ 

or 

►-*  c 

H 

►H 

cr 

» Q. 

• 03  ^ 

■K 

rvj 

V-  O 

<3^ 

cn 

lu  r:  lu 

OH  •'X 

• - — 

X 

★ 

■»—> 

X 

X' 

* 

LL 

in  + h* 

S3  Cr  r"  C 

'C  — 

m 

•k 

IL 

+ III 

in 

*-s 

1/5 

■ t — • in 

•— • ■<  • O. 

in  • 

y 

<»x» 

37  *x 

2: 

W 

rg  n 

LL! 

2- 

c c. 

HHC  7 

oj  cc  eg 

” 

CC  C 

— 

Q 

+ 

u in 

cc  < 

37  • H 

CC  10  * *-• 

x in  33 

in 

33  7* 

M 

+ c 

* 

M 

C: 

X ^ . (T 

•<  7".  ifl  *Z 

CM 

l— i 

x H 

M 

H 

3* 

• 

• 

HH 

w m rg  <; 
CL  'O  — _! 

f t • 

X + 

Lu 

Or 

U'  LX 

• 

X C 

— >' 

rg  rg 

u: 

*z*. 

C 

U’ 

in  ix  ^ uj 

* -)  O 

c 

a o 

mx^ 

cr 

7^ 

• * 

X X 

in  x 

>C 

r. 

cr 

in  ~ m 

» — < • 7-"  'O 

X 1-1 

x < 

IM 

a 

<* 

r\j 

•— « • 

M 

1 

LU 

in  >-  ^ cr 

cc  •»  m *-«  m 

— X 

— 

in 

“ 

pi  cr 

• 

hH 

— 

rn 

— 

pg 

PM 

C 

i— * 

< < u : 

IX  H * H rg 

- II 

• 

X.* 

• 1 1 ’ 

a 

X 

• 

XX  • 

»■ 

» 

* 

*/ 

• 

C: 

m oi  (— • 

H x eg  • 

*~>  ^ 

c* 

pfc- 

o •> 

in 

r> 

f-* 

X 

w H 

M 

»— « 

V/ 

r- 

in 

o 

~x 

37  & 2 37 

H 

77  < •Cr  — 

W 

III  cr 

*— « 

m < 

— 

n ri 

* » • 

H O 

N— ^ 

| 

X 

— 

• 

• 

11! 

!x  H — < II 

LL 

• 

H < 

• 

It 

+ 

■x 

• 

IX 

II 

CC  • 

X '1. 

C 

* n 

X** 

IL 

l1  ' 

ii 

o a 

• 

X 

H — s 

< 

= in  v-  _j  -> 

II  o H 

m 

H IX 

II 

H cr. 

M 

c o 

II 

II 

r° 

< 

II 

II 

i— * 

n 

ii 

t-4 

33  33  in  •• 

_J 

••37  33 

x in  in  or 

in  t-< 

• 

LU  (C 

in  '.: 

•»“% 

*-N 

7* 

*"X 

«*■% 

m u 

^ i. 

w w 

•« 

H ••  C O > 

— rvi  < U1 

< < c o to  in 

■-C  >-  O 

n x- 

r\i 

c* 

V.rf> 

X 

i— » 

c 

— 

rg 

in 

Ct  y.  UJ  S',  o 

3CH7- 

*-> 

• + 

_l 

M 

i 

j y 

M 

rv 

n 

*— ♦ 

fL 

P^ 

o >< 

» 

•• 

H 

x-» 

• 

m 

H 13 

» 

• 

c\ 

SJ  ?H 

•< 

—*  < <• 

II 

*->  “> 

x*p 

cr 

. — >— i 

IT 

c 

1-^ 

<1* 

hH 

II 

1— » 

H- 

*r 

1— 1 

x-» 

U' 

r r«au.crujc«c 

*-> 

W > 

IX 

T-' 

•** 

X 1- 

r) 

V. 

• ■* 

X 

r~* 

5 'JL 

N-X 

— 

x^ 

x-» 

X- 

x-» 

x-» 

? ! 

> 

'nnnaaacroticxiUL 

•x 

c 

o 

m n (/;  c 

'O  V 

rr 

c. 

X *— 

IL 

!l. 

O 

LL 

tL 

1 * 

13 

r x x 

•— « 

7" 

H 

1 

X 

! 

tf 

l 

| 

! 

1 

> 

1 

i 

i 

c 

j 

*. 

O 

c 

1 

i 

j • — 

. — 

1 

rj 

i 

r\j 

i 

v*- 

p) 

:n 

i 

(J 

O 

o 

o 

o 

C c:  n c o 

r.  o c C:  r- 

p 

G'  C. 

o 

r-) 

r: 

occ 

C 

^ * 

c’ 

* ' 

c 

f 

X 

C' 

r • 

v* 

(L 

f' 

fV> 

r“ 

— rg  rg  m rr 

uTO  r*  3 > 

c 

— eg  <n  m 

.n 

c r- 

33  O l 

— 

nj 

p^ 

c 

r- 

> 

C; 

— 

p;  O 

P' 

n 

O^OO 

o O O c~  c 

— 



— 

•— 

— _ 

_p 

— 

n* 

r»! 

n • 

r'  i 

r' ' • 

r • 

**’j 

r'  • 

n' 

•Y' 

r ■ 

r 

'■** 

rv> 

ro  (vn 

<'0 

*r> 

p^  r**)  < ^ 

'0 

rp 

V 

r** 

f'T- 

r 

(v 

P 

(v 

( 

r' 

P' 

r * 

B-6 


* 

* 

X > 

* 

■X 

X or 

* 

X 

X UJ 

* 

■X 

X f- 

* 

X 

X 

* 

■X 

X tO 

X 

X 

X :? 

X 

•x  r- 

X c 

♦ ' — ■ 

x — 

X _l 

* o 

x tn 

X 

X ~ 

-x  r 

X Cr 

x c. 

X O: 

X c 

x c 

X Li' 

X X 

* T. 

X (-• 

X 

* 1- 

X 

X H 

X Li  ' 

x r* 

X c 

* r= 

X ! J J 

X zxz 

* 

X Cr 

x H 

* f 

X X 

X 

* rz 

X 

X LU 

* r 

X X 

X V 

x ' 

X CO 

X *~i 

x H 

X W 

X _l 

* a. 

x zr 

X 

x r 

X 

X CO 

x 

X t-1 

x r 

* X 

X ZO 

X — 

* in 

X c 

X H 

4:  •— i 

X 

X zz 

* H 

X UJ 

X co 

* nr 

X t- 

X >— • 

x zo 

*-> 

X < 

X UJ 

* CJ 

<M 

X o 

OJ 

x ?; 

■X 

■X 

X 

X 

X 

* c 

■X 

X 10 

X 

X X 

* H 

Ul 

X ZZ 

LU 

X H 

■X 

on 

X _J 

to 

X •— 

* o 

l-H 

X CL 

M 

x v: 

-X  2 

'*■  /->  ^-x 

X 

c ^ ^ 

X 

-X  •— < 

L"  UJ 

X T 

5^.  U’  UJ 

X >c 

~ fC 

X O:  Cr 

X 

x cr  a 

X 

m 

* 

cr 

+ - * 

X 

+ - . 

X 

*r 1 

N-m> 

X 

n 

c\j  <m  t?  r 

X I-1 

CM  <M  X T 

X 

»— « 

^ ' 

X 

cj 

x (->  t—  r 

X 

X C 0 *"■* 

X 

cr. 

X 

o 

X-  < >'  X O: 

x m 

x < :x  x 

D* 

x 

•n  x 

•»_ 

X 

< 

N P”  V-*  V-*  ^ 

x < 

S'*  *-• 

v^- 

X 

< CJ 

x 

X 

CM  X CM  CM  X 

X 

(VXfNJfV 

X 

x 

LU 

in  *- 

X 

t:  zt.  to 

**■> 

X Z' 

#v/  7:  ^ 

(P 

X 

z-‘  X ro 

7"  — 

X 

•-  S.  < < ZZ  2-. 

s 

X r3 

s.  < < 

Z' 

•*T 

«■ 

s 

X 

oax 

G 

< c 

X 

n 

C\J  H*  fr— * G < C 

X Cr. 

^c\HH 

X 

< 

O 

X 

c in  X 

r* 

/"A 

X 

H 

xx  •■:  < z~  c 

— 

X H 

X X < < 

ZZ1 

— 

X 

1—  CM 

l—l 

H V 

X 

CJ 

+ X ^ H1  1— ■ 

x o 

+ X 

X H 

X 

OJ\ 

U! 

X 

U! 

CM  CM  00  ZO  L! ! 

X Li 

cm  n1  in  zz  cz 

LU 

J*. 

x 

m < • 

Cr 

Ct 

X 

Q. 

X O C *-i  Cr  O 

i-h 

X X 

X O c •— i <•• 

Cr 

a 

l-H 

X 

a.  r — 

N 

1"  — J 

X 

LO 

<M 

X < cj  m ~ x 

U 

X 00  CM 

X < cj  tn  m 

LL! 

U 

X 

00  or  li 

vt 

t-H  ” 

X 

— 

^ X X X 'C  *—* 

z 

X 

r - : xx* 

r~ 

l-H 

s 

x 

o 

• 

Cr  w 

X 

Cr 

P 

— x co  co  «n 

v-» 

o 

to 

— X CO  CO  rf 

• 

Cr 

x-*c 

X 

Cl  T~.  CM 

/■" 

7 - 7*. 

X 

f i • 

. — < < c zj 

X f: 

• 

* «c  <c  in 

o 

rr 

x 

•'  • 1 

C?  •-* 

X 

T” 

— 

^ r\» 

Hhrzuc 

l-H 

X 

— r\j 

•-I  t—  z z-  — 

• • 

X 

l-H 

X 

- 'L  < 

• 

*.1.  < 

X 

C* 

II 

— » 

ru  x x • u. 

< 

x r 

li  — • 

^ CZ  X X II 

c o 

• 

It, 

< 

X 

C"  c 

H 

2J 

X 

a. 

l-H 

u.  c ii  ii  f- 

IT 

X CL 

l-H  » l-H 

X CM  II  *— 

II  II 

H 

— 

X 

a.  • 

00 

t: 

r % 

V. 

> 

l-H 

>— » 

M 

00 

kA-s 

s 

CP 

i • 

C.J 

v’ 

l-H 

ii 

fp 

A.  /A 

00 

Ul 

CJ 

V/ 

x x o 

n r 
• \j 

m 

CL 

o 

1 L 

CV* 

II 

— 

p.» 

< 

•pi 

cr 

* 

X o 

LL 

ii 

— CM  C 

— CM 

<r  in 

cc 

X X CO  II 

» 

» 

-^v 

1 

7^ 

<- 

u. 

ii 

« > 

CO 

• 

• 

C’ 

X 

* 

1 1 

r-  < u 

n 

» A 
•X 

CO 

» » cr 

m • 

, j 

Cl 

t—' 

X X < a 

t—t 

’w' 

' 

u 

X — ' 

II 

<r" 

•r 

l-H 

l-H 

L'l 

X 

l-H 

ii 

- 1 

-c 

**r 

•-H  l-H 

•-H  *-H 

Ax’ 

u. 

X— 

j- 

X •_  Cc  — 

>*— * 

i- 

> 

■ ~ 

1 ' 

X O' 

I'J 

l-H 

7 

x-<* 

5-L 

> 

u. 

-V 

r. 

C 1 ” 

l-H 

^ . 

VX  A/  C 

w *—* 

X > 

• J 

I'J 

x n u ^ 

1 X 

i*  « 

■*  > 

X c 

7s 

a* 

X 

* 

’I- 

*-H 

2 

o 

■H 

JO 

n x- 

X 

‘v 

x ll  n x x 

l-H 

7- 

L? 

CJ 

x r > < 

l-H 

*J<  r 

►— 1 

> 

l-H 

< 

1 

X | 

1 

4? 

1 

| 

•V  | 

1 

X 1 

1 

-V 

1 

1 

X 1 

1 

X 1 

e 

1 

-i< 

1 

a 

C' 

1 

X 1 

;> 

1 

m 

X 1 

-c 

1 

> 

1 

r- 

O' 

CP 

1 

X 1 TO 

M 

CJ  CJ 

CJ 

CJ 

JJ 

CJ 

O CJ  CJ 

r 

" 

' s 

r 

C ' ' 

c 

r~-. 

c 

c 

' ‘ 

n » 

r 

— 

p 

(“ 

i*  • 

c" 

r rN  — 

. . .V 

f“ 

~ * 

" *• 

c 

' . O CO 

'•j 

,n 

r~ 

rp  O 

— 

r\j  rp 

M 

JO 

>o 

r~ 

r a 

O' 

— 

p^ 

n. 

•c 

r~  co  c>  n 

pA  > 

C. 

— 

CM 

p 

•j  o c o 

XT* 

rv»  nn 

-X 

*■  r 

» 

* * 

•t 

\ »- 

-*■ 

, 1 

A 

■*« 

A • . 

i "* 

A 

•A 

* o ^ 

A «> 

,“A 

A 

.A 

■“ 

r . 'P  p 

P 

f* 

r 

f 

AT 

P P 

CO  P^ 

P*»  rp 

P' 

r*' 

CC 

ro 

p' 

r-* 

xv* 

ro 

P.  P' 

P' 

p: 

P 

O ' P P 

p 

P' 

n> 

P « 

p 

O'  p p 

■nUi 


4 


-An.- 


'~n  f\a— 


~~r  7 


f\i 

* 

* 

fM 


fM 

* 

* 

:.u 

in 


+ 

fM 

* 

* 


a- 

+ /-**  <«. 
fM  ^ ^ 

* UJ  LU 

* a:  a: 


^*N 

fM  ♦-* 

T 

• • 

* 

**** 

M 

*— * 

y.-‘ 

♦ 

* 

< 

X 

X 

Q 

•— * 

v-» 

m 

fM  V 

r\j  cm 

X 

m 

a. 

~ 

7-.  71 

CO 

in 

— 

fM 

»— » 

fM 

< < Z 

2:  “ 

• 

rn 

* 

• 

H • 

in 

s-*- 

a 

h*  H Ct  < C 

n 

in 

* 

r- 

a-  m- 

\ 

LL  < 

< 

< m 

cr  — 

«■■**. 

• 

‘O 

r- 

c o 

ULI 

+ 

V— ' 

H * r~*  X 

m 

r- 

> 

O 

!M 

fM 

in 

CO  — 

in 

(M  X 

or 

7~  LL 

2* 

N-* 

in 

— 

ii: 

is- 

* 

\ 

fM 

+ v.  r-  —<  c 

fM 

* 

•K 

•—  O; 

C *-i 

m 

< 

m 

* 

or 

in 

* 

<^*N 

* 

iu  in  m 

c n 

* 

fM 

o m ^ 

UJ  — 1 

* 

• 

'U 

fM 

5> 

in  in 

OJ 

7"  fM 

►— » 

in 

♦ 

* X 

*— * z 

* 

• 

f\» 

on 

• 

\ 

• 

n 

m 

»— i » — < 

* 

O:  * 

’>*>■ 

* 

o cn  • 

O 

• 

fM 

• 

rv 

— 

• 

*+** 

co 

— • . 

n n 

— 

M fM  < 

< < o 

X'  z 

PJ 

— 

<c 

11 

1 

o 

li 

< 

• 

ii 

7*:  r* 

II 

^ II 

n 

ii 

f\J 

»— « 

»— 1 

5?" 

T LL! 

Cr  * — 1 

\ \ 

II 

M 

•— « 

c 

+ 

H O 

*— • 

Or  Cr 

»— * *—* 

<^*N 

♦— ♦ 

— 

• 

X 

X • 

LL  < 

• 

• 

* — ' 

• 

X 

in 

w - 

Cc 

n 

II  II 

*— < 

»— i 

LL 

< 

11 

II  H 

— 

•— 

o 

• 

o 

OliiC 

III  S’J 

c 

< o 

1 

c. 

II 

v-o 

^ CO 

LL'  6 

n 

n 

o 

n 

o 

II 

ii 

o 

— 

in 

m in  fM  in  in  m 

II  ^ 

M* 

•o 

LL  fM 

II 

f\i  <t: 

in 

O: 

*"*N 

*"*N 

o 

— 

— ■ O 

*"*N 

ii 

— 

ji 

11 

*— i 

— 

t«4  ►— » 

— 

— 

ex- 

II 

o 

• 

• 

C:  — ! 

rn 

•+ 

*— 1 

ii  ir 

»-H 

rs 

♦— « 

in 

it 

T~ 

< 

< 

•— 1 

t-^ 

U — 1 

ho 

>*** 

S-*' 

c. 

i—i  or 

;n 

# ■ - 

7~*  71 

c 

s_ 

U 

t—i 

JS 

r ■ 

N— ' 

^ IL 

> < t 

LU 

•-t 

< < 

ci 

<nx:i 

< in 

in 

c/>  c:  n a:  n: 

< cr 

< O O’ 

X 

V 

X 

'X. 

X •— 

z o 

I — l 

cr  mi 

o 

o 

C/ 

c 

o 

1 

1 

o 

r 

— 

<M 

ro 

tj- 

in 

1 

> 

** 

— • 

1 

o 

c 

c 

c o o r.- 

o 

c c 

c 

p, 

o 

c T 

r-'' 

c c 

c 

c 

r 

c: 

C' 

c; 

c r> 

•*  ' f 

•*> 

V 

x O' 

c 

— 

r\jfOTf 

in  'C  o- 

#v»- 

>A_ 

> 

C 

— 

fM  "n  m- 

in  o n- 

«v 

O'  c 

— 

CM 

n 

V lO 

:c  r* 

>> 

> 

'T 

•■r 

r- 

r*»  ^ 

r*- 

r^. 

r- 

nr> 

f>n 

<*•  rr 

or 

'r> 

N 

N 

\ 

\ v 

% 

> 

ro 

ro 

ro  ro  r« 

m 

rr 

<o 

n 

c* 

m 

>y- 

CO 

ro 

r*“ 

fO 

0* 

CO 

m 

/nr>  fV' 

rr 

po 

B-8 


„ --  ■*  ■.*.  — .».«■■■> 


CJ 

X 

« 

H 

m 

n 

ai 

4= 

Q 

§ 

0 

4J 

00 

c 

Cm 

<0 

0) 

HQ  <M 

H C — 
a.  *-i  in  — 


in 

_J  in 

O O C 

2 O 8c  8 

Sr 

M 

Q C 

U < COHO 

_J  cc 

OP  > 

» 

a c r~  o cc  rr 

a c 

<M  UJ 

nj  Q:  *-i 

• ccaaocpou. 

* H 

m cr 

1 H UJ 

OHUlHCCaOC 

C Ul 

c^>  <x 

II  in  O 

ll  in  o in  ti  ii  < it  m 

ll  a 

a)  § : 
o i-i  o 

MH 
3 0 = 
O In 

m a ^ 
ai 

10  x> 
•H  B 
x:  3 
H C 


u. 

uj 

no  uj  c c < < 

2vctO(/5tn(ni/iUi><<Lijnacro:in<<D.!ija.Lu<< 

j-ji/Koir'.oi/H/XQi-'rHUJii.'rQsnNrMQQ. 

wm  x>  x)ca.'vm=>in~j'S)~  a \>  n ?.  _j<ohuh_jh 


ll  u.  cr  cr 
n c c < -tj 

n uj  xi  o.  H 

ooccoomoooooocooooocooco 

-■Mi*'tinr,flN'T';\0-Mr»'iti(vcKT:vir  — -m  ■r> 
OOOOOOOOOO  — — — — — — — — — — <Mf\j<N<\l 


2-iOAbORT  LDA  =1 ,DL 

2o0  S1A  IPARTY 

260  IRA  RETURN* 1 

2 70E0F  LDA  = 1,DL 

280o5T'JRrJ  STA  I EOF 


listing  of  the  program  "WTB",  magnetic  tape  write  program. 


Source  listing  of  the  program  "WTBZ” , magnetic  tape  write  program. 
This  program  writes  a record  of  length  1 word  onto  tape  with  assigned 


Batchjob  listing  of  the  sequence  of  operations  to  be  performed  in 
copying  one  magnetic  tape  (assigned  name  "IN")  onto  another  magnetic  tape 
(assigned  name  "OT")  which  may  be  read  without  parity  errors. 
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9220C TEMP  TEST 

9230  SUBROUTINE  TESTI 

9240  COMMOn  /BUPFER2/  F(5 12.2 ) ,G( 51 2 .2) . A(5I 2 .2 ) 

92bO  CALL  n'rt 


0260  CALL  FFT 

92  70  CALL  IFF! 

9280  RETURN 

9290  END 

9300  SUBROUTINE  TEST2 
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9600C SUPPRESS  LOW  TEN  ELEMENT?  TO  E0UIL  TWICE  FORMANT  PEAK 

9610  COMMON  /BUFFER2/  F( 6 1 2 . 2 ) ,G< 5 1 2 , 2 ) . A (5 1 2,2 ) 

9620  FMAX=S0RT(F( I . I )**2+F( 1,2) **2 ) 

96 JO  DO  100  1=2,221 

9640  100  FMAX  = AMAX» ( FMAX ,SQRT(F( I, l)**2+F< I.2)**2  )) 


9650  FACf=SQRT( F( 257, I ) **2+F< 257, 2) **2 )/FMAX 

9660  00  110  1=2^3,266 

96/0  F( I, 1 )=F( I . I )*2/FACT 
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